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Resumen 


En el presente trabajo se exponen las herramientas utilizadas en en estudio de 
la estructura de bandas del grafeno y nanoribbons de grafeno con condiciones de 
borde zigzag (ZGNR), armchair (AGNR) y barbados (bearded) dentro del modelo del 
enlace fuerte (tight-binding). El transporte electrónico en nanodispositivos de grafeno se 
describe dentro del formalismo de Landauer y Landauer-Búttiker, dentro de las cuales 
usamos como herramienta principal las funciones de Green. Consideramos que cada 
nanodispositivo se encuentra formado por el nanodispositivo conectado a contactos 
semi-infinitos ideales, reduciendo en estos casos el efecto de los contactos mediante el 
uso de las autoenergías que se calculan usando las funciones de Green de los contactos. 
Los contactos de los nanodispositivos se consideran formados por cadenas verticales, 
donde cada cadena interactúa con las cadenas paralelas y vecinas, esta propiedad 
permite utilizar el método recursivo de convergencia rápida para el cálculo de las 
funciones de Green de los contactos y el método de Decimación para el cálculo de 
la función de Green de los nanodispositivos. Evaluamos la estructura de bandas y la 
conductancia eléctrica de los nanoribbons de borde zigzag, borde armchair y borde 
bearded, encontrando que los nanoribbons ZGNR tienen comportamiento conductor, 
los nanoribbons AGNR tienen comportamiento semimetálico si M = (N + 1)/3 (con 
N número de dímeros por cadena unidad) es un número entero y en caso contrario 
los nanoribbons AGNR tienen comportamiento semiconductor donde el gap de energía 
prohibida AF tiene relación inversa con el ancho del nanoribbons. Los nanoribbons de 
borde bearded tienen comportamiento conductor al igual que los nanoribbons ZGNR, 
en cambio los nanoribbons de borde mixto zigzag y bearded tienen comportamiento 
semiconductor, donde el gap de energía prohibida disminuye a medida que aumenta 
el ancho del nanoribbons. Analizamos también el transporte electrónico en nanoporos 
rectangulares en redes cuadradas y en redes de grafeno, calculando asi mismo el efecto 
sobre el transporte electrónico originadas por moléculas de diferentes configuraciones 
ubicadas en el centro de los nanoporos, los resultados muestran que las moléculas 


producen la presencia de resonancias y antiresonancias en diferentes valores de energía 


TI 


que son caracteristicos para cada tipo de molécula en el espectro de coeficientes de 
transmisión, estas resonancias y antiresonancias son evaluadas en términos de variación 


del coeficiente de transmisión. 


Palabras clave: Grafeno, tight-binding, transporte cuántico, transporte balistico, 


funciones de Green, formula de Landauer, nanoribbons, nanopores 


IV 


Abstract 


In this work we show the tools used to study the band structure of graphene 
and graphene nanoribbons with zigzag edge (ZGNR), armchair (AGNR) and bearded 
within the tight-binding model. The electronic transport in graphene nanodevices is 
described within the Landauer and Landauer -Biittiker formalism, within we use as 
main tool the Green's functions. We take in consideration that each nanodevice is 
formed by the nanodevice connected to ideal semi-infinite contacts, in such cases we 
reduce the effect of contact by using the selfenergy that are calculated using the 
contacts Green's functions. The contacts of nanodevices are formed by vertical chains, 
each chain interacts with the parallel and neighbors chains, this property allows to use 
the recursive method for quickly calculate the Green's functions of contacts and the 
decimation method for the calculation of the nanodevices Green's function. We evaluate 
the band structure and electric conductance of nanoribbons with zigzag, armchair and 
bearded edges, finding that the ZGNR nanoribbons have conductor behavior, AGNR 
nanoribbons have semimetallic behavior if M =(N + 1)/3 (with N number of dimers 
per chain unit) is an integer and otherwise the AGNR nanoribbons are semiconducting 
behavior where the energy gap AE are inversely related to the nanoribbons width. 
The bearded nanoribbons have conductor behavior as the ZGNR nanoribbons, however 
nanoribbons with mixed edge zigzag and bearded have semiconducting behavior, where 
the energy gap decreases as the width nanoribbons increases. We also study the 
electronic transport in rectangular nanopores on square networks and graphene, we 
calculate too the effect on the electron transport induced by different molecules located 
in the center of the nanopores, the results show that the presence of molecules produce 
resonances and antiresonances in different energy values that are characteristic for 
each molecule type in the transmissión profile, these resonances and antiresonances are 


evaluated in terms of vatiarion of the transmission coefficient. 


Keywords: Garphene, tight-binding, Quantum transport, Ballistic transport, 


Green's functions, Landauer formule, nanoribbons, nanopores 


Índice general 


1. Introducción 


2,1. Consideraciónes de TRÚ| + «39 2... 95:20. e a A 
2,2. Esimejura de Dadas vs ori ea e 
2.3. Aproximación en puntos de Dirac| ..................... 
2.4. Grafeno en un campo magnético uniforme] ................ 
CEA aaa 

2.5.1. Red cuadrada finita en campo magnético uniforme ....... 

O 
o ea 


261 Red cuadrada... 


APM 


3. Transporte electrónico 


3.1. Formalismo de Landauenh 


VI 


3.2. Formalismo de la función de Greedl ............. o. 29 


32.1, Función de Greeaal o 30 


3.2.2. Funciones de Green fuera de equilibrio] .............. 31 





oeficiente de transmisión y la función de G 34 
3.2.4. Método de decimación 36 
3.2.5. Funciones de Green de contacto derecho! ............. 39 


3.3. Formalismo de Buttiker para 4 contactos] . ................ 42 


4. Nanoribbons de grafeno 48 


4.1. Estructura electrónica de nanoribbons de grafeno zigzagl. ........ 50 
4. T.T. Método del Hahtbndidals osos sarr naa 52 
4,2, Estructura electrónica de nanoribbons Armchait| ............... 54 
4.21. Método del tight-bindinal. ¿sic nonas a 57 
4.3. Nanoribbons barbados (bearded)| ...................... 59 


4.4. Nanoribbons zizag-bearded|. .......... o... .. . . .. ... 0. 61 


4.5. Conductancia eléctrica... 62 


4.6. Nanoribbons en campo MagnétiCO| ......... o... e... o... 68 


75 
5.1. Nanoporos en red cuadrada] ............. o... ....... dí 


BILL Moadela. 78 


a.1.2, Resultados Y QISCUSION) e. 64 sie mo e me e BR 82 


VII 


6.2. Nanñoporos de pralenól. ¿2 «cien css e a 


BALL Elmoadea 


5.2.2. Resultados y discusión 


6. Conclusiones 


. Bibliografía 


A. Relaciones importantes 


A.1. Elementos de matriz para el hamiltoniano del grafeno| . .......... 


A.2. Relación de dispersión |f(k) 


VII 


101 


104 


112 


Capítulo 1 


Introducción 


En el desarrollo de la industria de los semiconductores de silicio, los circuitos 
integrados consisten de transistores, que hacen que el estudio de la electrónica y 
microelectrónica tenga importancia cada vez mayor, áreas en el que el transporte de 
electrones se describe por características estadísticas promedias, siendo las cantidades 
físicas utilizadas cantidades macroscópicas, ello debido a que el tamaño de los 
dispositivos es mucho mayor que la longitud de coherencia de fase de los estados 
electrónicos. De acuerdo a la ley de Moore el número de transistores por unidad de 
área en un circuito integrado se duplica aproximadamente cada dos años, el tamaño 
de los circuitos integrados caerá a 30 nm en el siglo XXI, tamaño en el que los 
circuitos integrados ya no pueden ser descritos por leyes macroscópicas, debido a que 
las características ondulatorias y los efectos cuánticos cobran mayor importancia, hecho 


que justifica el desarrollo de la nanoelectrónica. 


En la naturaleza encontramos que los átomos de carbono constituyen un elemento 
esencial de la vida en el planeta y la química orgánica, se encuentran formando 
diferentes compuestos de diferentes estructuras y diferentes propiedades. Así mismo 
el carbono puro en la naturaleza se encuentra en dos formas: grafito y diamante, 
que tienen propiedades eléctricas y mecánicas diferentes debido a que sus estructuras 
cristalinas son diferentes. Mientras el diamante constituido por enlaces híbridos sp*, 
es un material aislante muy duro; el grafito formado por enlaces híbridos sp?, es 
un material conductor y a su vez lubricante. Los primeros alótropos de carbono 
encontrados artificialmente fueron los fulerenos en 1985, que son moléculas esféricas de 


átomos de carbono considerados como sistemas cerodimensionales, este descubrimiento 


originó una intensa investigación de posibles materiales puros de carbono. En 1991 se 
realiza el descubrimiento oficial de los nanotubos de carbono con paredes multiples 
MWCNT y en 1993 nanotubos de carbono monocapa SWCNT, los nanotubos son 
sistemas unidimensionales, debido a que su relación longitud - radio es muy grande. 
En el 2004, en la Universidad de Manchester, el grupo liderado por Andrei Geim y 
Konstantin Novoselov obtuvo experimentalmente, por primera vez, una monocapa de 
átomos de carbono sobre la superficie de un substrato de dióxido de silicio SiO» (1112), 
a esta monocapa se le denomina grafeno. Una sucesión de grafeno forma el grafito, 
el material más conocido por el hombre desde que se encontraba en las cavernas. La 
técnica desarrollada por el grupo de Geim es relativamente simple y está considerada 
dentro de las técnicas de escisión micromecánica, y consiste en pegar una cinta adhesiva 
sobre un pedazo de grafito de óptima calidad para luego colocarla sobre un substrato 
de dióxido de silicio. Luego de una búsqueda con un microscopio óptico y electrónico 


se detecta el grafeno en monocapas, bicapas y terrazas de varios niveles. 





Figura 1: Andre Geim (iquierda) y Konstantin Novoselov (derecha) ganadores del Premio 


Nobel de Física 2010, por el descubrimiento y caracterización del grafeno. 





Figura 2: Construcción teórica de otros alótropos del carbono a partir del grafeno: 
fullerenos (0D), nanotubos (1D) y grafito (3D) [1). 


Durante muchos años se presumía que sistemas bidimensionales puros no se 
encontraban en estado libre en la naturaleza. Opinión que se fundamenta en la 
argumentación de Landau y Pierls, quienes demostraron que cristales bidimensionales 
no pueden existir, por ser estas termodinámicamente inestables [8H6), dicha teoría se 
puso en tela de juicio con el descubrimiento del grafeno en el 2004. Sin embargo se 
puede considerar que las teorías de Landau y Pierls siguen siendo válidas si se tiene en 
cuenta que las hojas de grafeno presentan ligeras corrugaciones y defectos en la tercera 


dimensión y por tanto no son estrictamente bidimensionales [6117]. 


En el grafeno, los átomos de carbono se unen químicamente formando una 
estructura cristalina bidimensional de tipo hexagonal, que resulta de la hibridación 
de los orbitales 2s, con los orbitales 2p, y 2p,, que produce tres orbitales híbridos sp? 
que se encuentran el en plano xy formando entre sí 120”. La hibridación sp? proporciona 


la estabilidad de la estructura bidimensional y las propiedades mecánicas del grafeno. 


De acuerdo a la propuesta de L. Pauling [8| los tres estados híbridos sp? son: 








mo >. v2 
25Pa) = B 25) + B 2px) (1.1) 
A 1 
asp) = 125) — —/2p2) — —=l2p4) (13) 
Pe v3 v6 Pe 2 Y l 





El electrón en el orbital p, perpendicular al plano de grafeno permance sin hibridar y 
le permite al grafeno adherirse al susbtrato de Si03 mediante fuerzas de Van Der Waals y 
simultáneamente proporciona al grafeno unas propiedades electrónicas singulares que 
ha llamado la atención de la comunidad académica mundial y que ha generado una 


intensa actividad de investigación tanto teórica como experimental. 


La estructura electrónica del grafeno dada por los orbitales p, presenta dos puntos 
K y K”, en la primera zona de Brillouin, donde las bandas de valencia y conducción 
coinciden puntualmente, los portadores de carga en dichos puntos no poseen masa y la 
relación de dispersión de la energía en las proximidades de los puntos K y K” es lineal 
respecto a la magnitud de vector de onda k, es decir E = +hkvup, donde vp = 10% 
m/s, es la velocidad de Fermi, una velocidad que es 0.33% respecto a la velocidad de 
la luz, siendo considerado casi relativístico, esta propiedad es una de las razones por 
el cual se considera al grafeno como el primer material de estado sólido para realizar 
experimentos de electrodinámica cuántica, que permite probar por ejemplo, la paradoja 
de Klein. A los puntos K y K'” se les denomina puntos de Dirac y a los portadores que 


se mueven al rededor de dichos puntos se les denomina fermiones de Dirac. 


El comportamiento electrónico único del grafeno se debe a las características 
relativistas cuánticas inusuales de los fermiones de Dirac, estos bajo la presencia de 
un campo magnético uniforme dan lugar a un espectro de energía particular para 
cada nivel de Landau |9//12|, permiten también la observación del efecto Hall cuántico 
anómalo a bajas temperaturas y a temperatura ambiente [16]. 


Las redes finitas de grafeno, denominados puntos cuánticos, se caracterizan por 
el confinamiento de los electrones y los huecos debido al tamaño del sistema. Las 
colinas, ondulaciones y corrugaciones del grafeno afectan intensamente las propiedades 
electrónicas de los puntos cuánticos por ello es necesario poseer una lámina de grafeno 


lo más plana posible cuya obtención ya es posible [17]. 


En el capitulo 2 describimos la estructura electrónica del grafeno, la aproximación 
en las proximidades de los puntos de Dirac para evaluar el efecto del campo magnético 
uniforme sobre el grafeno. Se hace tambien una revisión de las mariposas de Hofstadter 


para redes finitas e infinitas en redes cuadradas y el grafeno. 


En el capítulo 3 desarrollamos el modelo de Landauer para el transporte electrónico 
a bajas temperaturas, describimos también la teoría de las funciones de Green para 
fines de calcular la densidad de estados y la conductancia eléctrica de nanodispositivos, 
evaluamos métodos numéricos que nos permiten calcular recursivamente las funciones 
de Green del dispositivo cuántico y de los contactos semi-infinitos, estos son el método 
de decimación y el método recursivo de convergencia rápida. Finalmente analizamos 
el formalismo de Landauer-Buttiker para evaluar numéricamente la resistencia Hall de 


un punto cuántico con cuatro contactos. 


En el capítulo 4 aplicamos la ecuación de Dirac para describir los estados de borde 
en nanoribbons zigzag y armchair, el método del tight-binding para la estructura de 
bandas de nanoribbons y aplicamos la fórmula de landauer para evaluar la conductancia 
eléctrica de los nanoribbons con diferentes bordes. Así mismo analizamos la estructura 


de bandas de nanoribbons bajo la acción del campo magnético. 


En el capítulo 5 aplicamos la formulación de Landauer al estudio de nanoporos 
de estado sólido construidos en redes cuadradas y el grafeno. Evaluamos tambien los 
efectos sobre la conductancia eléctrica producidas por la presencia de moléculas de 


diferentes configuraciones en el centro de los nanoporos. 


Finalmente presentamos el resumen de las principales conclusiones del trabajo. 


Capítulo 2 


Grafeno 


2.1. Consideraciones de red 


El grafeno es una monocapa de átomos de carbono distribuidos en una red cristalina 
bidimensional de tipo hexagonal, conocido también como panal de abejas. Cada átomo 
de carbono tiene seis electrones ocupando los orbitales 1s?, 2s? y 2p?. Los electrones 
1s? se encuentran fuertemente ligadas al núcleo atómico debido a su baja energía, por 
tanto no participan en los enlaces químicos. De los cuatro electrones restantes en los 
orbitales 2s? y 2p? de mayor energía tres producen enlaces covalentes hibridizados sp?, 
que proporcionan la estabilidad de la estructura hexagonal y hacen que los átomos de 
carbono se encuentren fuertemente empaquetados [1)[1sH20]. 


Los electrones en los orbitales perpendiculares al plano del grafeno son relevantes 
para las propiedades de transporte eléctrico, por lo que consideramos solo estas en el 
cálculo de la estructura de bandas de energía. La estructura hexagonal no constituye 
una red de Bravais, pero se puede describir por medio de dos redes de Bravais 
triangulares formada por átomos A (” e”) y B (” o”) respectivamente. De acuerdo a 


la figura [3] los vectores primitivos de cada red triangular del grafeno son dados como: 


3,.,3 
d, =vV3a4, di = as + 3% (2.1) 
siendo |d,| = |d2| = V3a la constante de red y a = 1.42 Á la distancia entre dos átomos 


de carbono próximos. 





Figura 3: Red cristalina del grafeno, representado por dos redes cristalinas triangulares de 


átomos A (” e”) y B (” 0”), donde a, y da son los vectores primitivos de la red real. 


Cada átomo de la red A en la posición R? tiene como vecinos próximos a tres 
átomos de la red B ubicados en las posiciones RF = RA + Ó,, j=1,2,3. De acuerdo a 
la figura [8] los vectores Ó, son las posiciones de los átomos B respecto a las posiciones 
de los átomos A, que se escriben como 


2 v3.1. 2 1, 


1 si 
01 = ———al — af), 0, = qe = 7%, 03 = 0. (2.2) 


Así mismo cada átomo de la red B en la posición RF tiene tres átomos de la red A en 


las posiciones Rf = Rf? — 0 como vecinos próximos. 


En el espacio recíproco, la red recíproca asociada al grafeno es otra red hexagonal, 


donde vectores primitivos b, y bz asociados a a, y dz son dados por 








2.2. Estructura de bandas 


Nuestro propósito es describir el transporte electrónico, por tanto nos centramos en 
esta parte a la estructura de bandas asociada a los orbitales perpendiculares al plano 
del grafeno, para tal fin consideramos que los electrones en los orbitales perpendiculares 
de los átomos en las sub redes triangulares A y B son representados por los estados 
cuánticos 19) y [62) respectivamente. La interacción entre vecinos próximos ocurre 


de modo que un electrón del átomo A (o B) pueda saltar a los tres vecinos próximos 


y 


B (o A), esta interacción se describe dentro del modelo del tight-binding, modelo en el 
cual el hamiltoniano tight-binding en primera cuantización se escribe como 


H ==) (1640?|+|07)(0?)) (2.3) 
(5,3) 


siendo t = 2.8 eV el parámetro de Hopping y (¿,j) una notación que establece que la 


suma en (2.3) se efectúa solamente sobre los vecinos próximos. 


Las funciones de onda que describen a los electrones en la red triangular A/B es 


dado por 
(162) = p(r — REP). (2.4) 


Siguiendo la propuesta de Wallace 24), las funciones de onda de Bloch Y y y? 
asociadas a las redes triangulares A y B se expresan en términos de las funciones 


de onda de los orbitales atómicos como: 


2 ARS (y — R?) (2.5) 
Da eri gr — RP) (2.6) 


teniendo en cuenta (2.4) reescribimos y como 

NA los) = (Y A) ) = (rl) 
j 

DA ros) = (E eE Os) ) = (eli) 
J 


de donde los estados cuánticos |W*) y |4P) asociados a las redes triangulares A y B 


ver) 


Y (1) 


tienen la forma 


pz) = RRA, (2.7) 


Jj 
Rp) = MEROS). (2.8) 
J 


Considerando (2.7) y (2.8) escribimos el estado cuántico |) del grafeno como una 
combinación lineal de los estados |1;*) y |yP), 


|eb) = cel) + co l0;). (2.9) 


Para calcular la estructura de bandas hacemos que el estado |) satisfaga la ecuación 


de Schródinger independiente del tiempo 


HlYr) = En Vx) (2.10) 


con el fin de resolver (2.10) expresamos en forma matricial el hamiltoniano A en la 


base (|), (194) 


(01H od) (we 1H (0) 
los elementos de matriz de A dada por la ecuación (2.11) se deducen detalladamente 
en el apéndice[A.1] y resultan 


(6H Ioé) = (11H 1w) =0. (2.12) 


mn deta a) Ca 


(1H 1uE) = 1109 2.13) 
donde f(k) = cikór sh cik-ó2 + ció 
Así mismo se encuentra 
id IO (2.14) 


Sustituyendo los resultados (2.12), y (2.14) en la ecuación (2.11), la expresión 


matricial de H es 
o f(k) 
H=- . 2.15 
e 0 ) e 


Los eigenvalores de energía de H dado por (2.15) se determinan de 


det(El — H)=0 


E tf(k) 
tf (k) E 





E" — 4 f(x) =0 
Ej = +t|f(k)| (2.16) 


siendo |f(k)| = /|f(k)1? = y/f*(k)f(k), obtenida en el apéndice y expresado 


como 
|F(«)] = | 1+ 4 cos? (Pra) + 4cos (Pra) COS (Gta) (2.17) 
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es decir la estructura de bandas del grafeno (2.16) es 





Ex = +Ht|f(k)| = +£,|1 +4c0s (Pe) COS (Gta) + 4 cos? (Pe) (2.18) 






¡de valencia 


Figura 4: Bandas de valencia y conducción del grafeno dada por la Ec. (2.18) 


La banda de energía Ej = —t|f(k)| representa la banda de valencia y Ej = t|f(k)| 
representa la banda de conducción. Cuando |f(k)| = 0 la banda de valencia y 
conducción coinciden, dichos puntos en el espacio recíproco se conocen como los puntos 
de Dirac K y K”. 


2.3. Aproximación en puntos de Dirac 


En esta sección verificamos que las bandas de valencia y conducción varían 
linealmente con el vector de onda en las proximidades de los puntos de Dirac K y 


K”, puntos que son dados de acuerdo a ([2.18) como: 
4 


K=—£,  K=-— T 2.19 
34 3a 3v3a ) 





expresando el vector de onda k como k = K + q, donde q = (q.,q,) es el vector de 


onda respecto a los puntos de Dirac podemos escribir f(k) como 


F(k) = F(K-= aq) a Aé — ya ¿óiaós (2.20) 


utilizando la aproximación lineal e? =1+=x 


só =1+iq:0, (2.21) 


A (1 iq: 55) =$ O ES 
si tenemos en cuenta que f(K) = 0 entonces 


F(k) = F(K) + SY ¿Kóiig . Ó; pt ¡Y eg 








ON aa 
F(k) =1 [> (434 = dy) + Fl > [V34s = dy) + aq,| 


donde e+i% = cos (2) +isin (2 1) = (- 1 +iv3) 


109 =35 [5 (1343) (V%0: 49) + 3 (14143) (V3no 01) +240 





su 
0) =i7 [V34s + qy + 1342 +1W3qy — W3Q% + Qy + 13q0 — 11 3qy + 4, 


Mesas, a ' 
F(k) = iz [6iq,, + 6q,] = 7 [—6q, + 169, 


3 
Ha) = ¿alg — 149) (2.2) 
con este resultado el elemento de matriz (4f|H (1%) se expresa como 
A B 3 . 3 ! 
(ÉIH JUE) =-+409 =+ ( Za(a: —i94)) =3atla.— ig) (223) 


del mismo modo el elemento de matriz (4P|H |) es 
._9 : 
(ALE ROE) = (WEA 14)" = zat(q, + ig) (2.24) 
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luego el hamiltoniano (2.11) en las proximidades del punto de Dirac K es 


3 0 == 
H = at .—00U 
2 Qu + Uy 0 





0 y 
H= hp OS (2.25) 
qe + iy 0 
3ta  3(2.8eV)(1.42 x 10m) á : E 
dond = = = 1 l : 
onde Up = 7 658 x LO160V-S 0” m/s es la velocidad de Fermi 


Resolviendo el problema de eigenvalores para el hamiltoniano (2.25) encontramos 





E —_ do 
det(E—H) = ci a 
—hvp (Qe + 1qy) E 
E? — Hui qu + qy)(4x — qu) =0 
E? — Puplaz +43) =0 





La Ec. (2.26) prueba que efectivamente en las proximidades del punto de Dirac K las 
bandas de energía E varían linealmente con el vector de onda |q|. Un procedimiento 
similar al rededor del punto K” prueba también que las bandas de energía en las 


proximidades de K' varían linealmente con el vector de onda. 


Los eigen estados correspondientes a la banda de energía E, = huplq| dado por 


(2.26), se expresan como |) = calW) + cal?) donde los coeficientes Ca y Cp se 


0 == 
Uz ar 14y CB CB 


(dr — 1Qy)cg =|lalca (92 +1qy)ca = [qlez 


de dy qe 1) E 
(e ie) eo=0s (Es + if) ca=cs 


(cos0 —isin0) cg =cA [(cosO +1sin0)ca = Cp 


obtienen de 


ca=e%cA 
asumiendo la condición de normalización |ca|? + leg]? = 1 y ca = o DS el eigen 
estado correspondiente a E, se escriben como 
1 


+) = 0 [Pé + 22H) (2.27) 
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donde tan Ó = y 
o 


De manera similar se encuentra que el eigen estado |, ) asociado a la banda de 


energía E_ = —hvup|q| se expresa como 


- A e710/2 — 619/2118 
|.) = 7! wi) ey]. (2.28) 


En el punto de Dirac K” se encuentra que f(q/) es 


/ 3 / el: 3 / 0d 
(q) = Uds - 1q,) = ¿Ud + 1q,,) (2.29) 
con lo cual las bandas de energía resultan 


El, = +hupla"| (2.30) 








siendo los eigen estados |y;,)asociados a las bandas E! dados por 








we) == lle Pty + e AP )]. (2.31) 


fl 
y2 
Considerando las matrices de Pauli o = (0,,, 0,) reescribimos la Ec.([2.25)) al rededor 


de K como 
Hx(q) = hur (079x + 049y) = RuFO - q = UFO : p. (2.32) 


De manera equivalente en la vecindad del punto de Dirac K” se encuentra Hx(q) = 


—Hl(q), de modo que el hamiltoniano total en los puntos K y K' es 


0 
== e oB ) (2.33) 
0  —o0*-p 


2.4. Grafeno en un campo magnético uniforme 


La dinámica de los electrones dentro de un material semiconductor convencional 
se describe usualmente utilizando la mecanica cuántica no relativista, dentro del cual 
se encuentra que la estructura de bandas de energía tiene una relación de dispersión 


parabólica con el vector de onda, de la forma E = Ey + h?k?/2m*, siendo m* la masa 





efectiva, Ey la energía mínima en la banda de conducción y E_ la energía máxima en 
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la banda de valencia. Bajo la presencia de un campo magnético uniforme de magnitud 
B la energía para el movimiento perpendicular al campo es cuantizado en un conjunto 
de niveles de Landau E, = Ex + hwe(n + 1/2), con n = 0,1,2,... y We = eB/m* la 


frecuencia ciclotrónica [25]. Este modelo falla cuando es aplicado a los electrones del 





grafeno, por tanto es necesario utilizar las ecuaciones relativistas de Dirac, considerando 
que los electrones del grafeno son fermiones de Dirac, caso en el cual de acuerdo a 
la estructura de bandas es proporcional al vector de onda, para analizar el 
comportamiento de los electrones del grafeno bajo la acción de un campo magnético 
B uniforme, asumimos que B = By, con By constante y teniendo en cuenta que 2 es 
un vector unitario perpendicular al plano del grafeno. Utilizando el potencial vector A 
tenemos que B = Vx A donde A = — Boyí es el gauge de Landau. La sustitución de 
Pierls [261127], Da —> P0x = Px — GA = Pr — eBoy, con q = —e carga del electrón, permite 
escribir el hamiltoniano como: 


y y Ú —P0w + Uy 
0 0 Doa WPy 0 
luego la ecuación de Dirac para H es 
0 Po=ipy 0 0 ve Ue 
0 +1 0 0 0 e y 
di 0 A ET ca CES 
0 0 0 —Pox + UPy Wiz Vic 
0 O —Ppow—ipy 0 vw, ye, 


donde 4% y 42 son los espinores asociados a los átomos A y B en K y Y£, y yÉ, los 


espinores de los átomos A y B en K. 


Las soluciones de 4% y 4% al rededor del punto K se obtienen de las dos primeras 


ecuaciones de (2.35), 
evt = UPÍ(DOs = py) VE 
EVr = UP(Pos + py) 
desacoplando 44 y y% de las ecuaciones anteriores, tenemos 
2% = Up (Pos — Py) (Pos + ipy)V (2.36) 
Vx = UilPos + Py) (Dos — ¡Py) VR (2.37) 
con la finalidad de encontrar los niveles de energía e para los estados 4%, consideramos 


que el hamiltoniano (2.34) es independiente de zx, es decir [H, p,] = 0, en consecuencia 
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p, es una integral de movimiento, [p., py] = 0, [y, py) = ¿h. 


E VR = UPlPho — 1PoxPy + ÍPyDos + PIJVR 
EVE = Up 
EVE = Up 


se 

e B 2 2 Pz á 21|,.,B 

(= E helio) Vk = [Bo (y a ) +Dy | Yx 
eBo 

Ef e Be E A 

=— |— 14heB Bo, 0 _ Pz El B 

2m (5 S .) VE 2m (y E) my 


Ve = (Zo yo) + E) k (2.38) 


(P» — eBoy)? — i[p» — eBoy, Py] + PER 
(p» — eBoy)? — heBo + pr) VR 











2 p2 
donde k= 2, y = E. 


El lado derecho de la ecuación anterior corresponde al hamiltoniano del oscilador 


armónico simple, oscilando al rededor de yy con eigenvalores hw (n + 3), con w = 


e2B? ] 
YE = 4/22 — £B0 es decir 
m m m 


1 E? eBo 1 
E A Ap, O 
Sm (q edo) 1 (+3) 





a 
= + heBy = heBy(2n) + heBo 
P 
e? = 2vñheBpn 
e = +up y 2heBon (2.39) 
donde n =0,1,2,.... como e tiene raíces positivos y negativos, podemos extender los 


valores de n a los enteros negativos y escribir los niveles de energía como 


e = sgn(n)y/ py In] = hw?“ sgn(n) y Inl, (2.40) 


: 2eB 
con wP tac = Ups] A 


Los niveles de Landau n pueden ser positivos o negativos, los valores positivos 





corresponden a electrones (banda de conducción) mientras los valores negativos 
corresponden a huecos (banda de valencia). Ademas dichos niveles no son equidistantes 
como en el caso convencional de materiales 2D, teniendo la máxima separación de 


energía AEmay = hw?" entren=0yn=l. 
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2.5. Mariposa de Hofstadter 


En 1976, D. Hofstadter predijo un raro efecto cuántico, conocido como mariposa 
de Hofstadter. El fenómeno consiste en que los electrones de una red bidimensional 


sometido a fuertes campos magnéticos muestran un espectro de energía complejo con 


características fractales [281/29]. 


2.5.1. Red cuadrada finita en campo magnético uniforme 


Consideramos primero el caso de una red cuadrada bidimensional finita de 3 x 3 
átomos, bajo la acción de un campo magnético uniforme B = By2, que es representado 
por el gauge de Landau A = —Bpy%, de modo que la sustitución de Pierls permite 
establecer que el parámetro de Hopping entre vecinos proximos sobre el eje x es dado 
por t, —> te*7M4, con m un número entero, y el parámetro de Hopping entre vecinos 
proximos en el eje y es ty —> t. En la figura [5]se observa la red finita de 3x 3 átomos, en 
la que el origen de coordenadas (x, y) se encuentra en la posición del átomo 1, de modo 
que el número de flechas paralelas al eje horizontal representan los valores enteros de 


m, que para el caso son m = 0,1,2. 


Figura 5: Red cuadrada finita 3 x 3 bajo la acción de un magnético uniforme. 


Teniendo en cuenta la sustitución de Pierls y el modelo del tight-binding el 
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hamiltoniano correspondiente a red finita de la figura [5] es 


0 t 0 t 0 0 ú uu 0 
t 0 t O. ¿ene 0 0 0 0 
0 t 0 0 0 ene yu 0 
ó 0 0 0 t 0 FU 0 
He |0 14" 0 t 0 t 0. ¿ete 1 (2.41) 
0 0 fee 1 t 0 UÚ UD 5 
0 0 0 t 0 0 0 t 0 
0 0 0 0 0 t 0 t 
0 0 0 0 0 E E 0 


El espectro de energía de que constituye la mariposa de Hofstadter 
correspondiente se muestra en la figura [6] donde a = 9/0, = p/q la razón del 
flujo magnético al flujo cuántico, es decir $ es el flujo magnético por celda unitaria y 
D. = h/e es el flujo cuántico, la razon p/q indica que a tiene valores racionales [29]. 
Siguiendo la misma idea utilizada para una red finita de la figura |5] calculamos el 
espectro de energías para redes finitas de 15 x 15 átomos y 25 x 25 átomos, cuyos 


resultados gráficos se presentan en las figuras [7] y [S] respectivamente. 


Energía / t 











Figura 6: Espectro de energía de una red cuadrada finita de 3 x 3 átomos en función de a. 











Energía / t 
o 





























Figura 7: Espectro de energía para una red cuadrada finita de 15 x 15 átomos en función 
de a. 
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Figura 8: Espectro de energía para una red cuadrada finita de 25 x 25 átomos en función 
de a. 
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2.5.2. Red ladrillo finito en campo magnético uniforme 


Con la finalidad de construir las mariposas de Hofstadter para redes finitas de 
grafeno, utilizamos la equivalencia topológica con la red de ladrillo, consideramos 
primero el caso de una red ladrillo de 3 x 3 átomos mostrados en la parte derecha de a 
figura|9] teniendo en cuenta que el campo magnético uniforme B = B/¿2 es representado 
por el gauge de Landau A = — Boy2 y que el principio de sustitución de Pierls 


hace que el parámetro de Hopping entre vecinos proximos es t, —> te" y by > 1. 





Figura 9: (izquierda) Red cuadrada de 3 x 3 átomos (derecha) Red de ladrillo de 3 x 3 


átomos 


El hamiltoniano tight-binding correspondiente a la red de ladrillo con 3 x 3 átomos 


representado en la figura [9] es 


0 t 0. tema q 0. 0. 0-0 
ii Dos 0 0 0. 0. 0-0 
0. t 0 0 O telir 9 0 0 
tera 9 0 0 t E 
H=| 0 0 0 t 0 to 0 terra q (2.42) 
0D 0 tetra 9 t 0. 0. 0-0 
0. 0.0 t 0 0. 0. t 0 
0. 0.0 DO tera 0 ft 0 4 
0. 0.0 0 0 0. 0. t 0 


Es espectro de energía asociado al hamiltoniano (2.42), que constituye la mariposa 
de Hofstadter se muestra en la figura [10] Del mismo modo las mariposas de Hofstadter 
asociados a redes finitas de grafeno constituidos por 15 x 15 átomos y 25 x 25 átomos 
se presentan en las figuras y respectivamente. 
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Figura 10: Espectro de energía de la red ladrillo de 3 x 3 átomos dado por el hamiltoniano 


en función de Q.. 





Energía / t 
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Figura 11: Espectro de energía de la red ladrillo de 15 x 15 átomos en función de a. 
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Figura 12: Espectro de energía de la red ladrillo de 25 x 25 átomos en función de a. 


2.6. Mariposa de Hofstadter para redes infinitas 


En esta sección reconstruimos el espectro de energía de redes infinitas bajo la acción 
de un campo magnético uniforme, que para el caso de redes cuadradas bidimensionales 
fue obtenido primero por Hofstadter ps), y verificado en varias publicaciones 36|, 


para el caso del grafeno 


2.6.1. Red cuadrada 


Sea el campo magnético uniforme B = B¿2, representado por potencial vectorial 
A = —Boy? de modo que la sustitución de Pierls introduce el factor de fase en el 
hopping entre vecinos próximos en el eje x dado por t, —> te?" y el hopping entre 
vecinos proximos en el eje y es t, —> t. La ecuación de Schódinger aplicada sobre el 


sitio (n, m) de la figura [13] produce 
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(n,m-1) 





Figura 13: Parte de una red cuadrada infinita en el que se denotan la posición del átomo 


(n,m) y sus vecinos proximos 


om == ud Obi CET =- re ad y MBE E ami en Unmadl (2.43) 
con Yr,m función de onda de los electrones del átomo en el sitio (n,m) y a = P/Do. 
con el fin de transformar (2.43) a un problema unidimensional consideramos 
ko = Ry =U0 Y 
Pr m= gr E Umm 


EUm == RN ES o + m1 03 tUm+1 


Em = 2t cos(27M0) Um + tUm-1 + tUm+1 (2.44) 


cuando q es un número racional p/q, p y q son primos entre sí, el tamaño de la celda 
unidad magnética es (qa, a), llamemos y(x+qa, y) = v(z, y) y v(x,y+a) = y(z, y) E 


2t cos(27ra) t 0 0 l t 
t 2t cos(4Ta) t | Ea 0 
0 t 2tcos(6ra) ... 0 dad 0 
0 0 0 “=. t 2tcos(27(q — 1)a) t 
t 0 0 “es Y t 2t cos(2179a) 
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Energía / t 





Figura 14: Mariposa de Hofstadter red cuadrada infinita 


2.6.2. Grafeno 


El espectro de energía del grafeno reproducido en esta sección, fué calculado primero 
por Rammal [87], el mismo que tiene la misma estructura recursiva que la mariposa de 
Hofstdater [28]. En nuestra deducción consideramos que el campo magnético uniforme 
es B = Bp2, que se representa por potencial vectorial A = —Bpy% de modo que la 
sustitución de Pierls introduce un factor de fase en el parámetro de hopping t.,. entre 
vecinos próximos en el eje x dado por t,, —> te" y mantiene invariable el parámetro 
de hopping entre vecinos próximos en el eje y, es decir £, > t. 
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Figura 15: Parte de una red ladrillo infinita en el que indica la posición del sitio (n,m) y 


sus vecinos próximos 


Aplicando la ecuación de Schródinger sobre los electrones de los átomos de las sub 
redes A y B situados en la posición (n,m) de la figura [15| tenemos 


EVn,m => ie E e lb. OA + A (2.45) 
En E tn dj RA a + ed ie E AE (2.46) 





con Unam Y Y; funciones de onda de los electrones en las sub redes A y B ubicados 


en los sitios (n,m) respectivamente y a = 9/00. 


Para transformar (2.45) y (2.46) a un problema unidimensional consideramos 
kz = ky =0 y 


Pam a prin => rn, Dr E e Hana ik mat - Y 


EUm => ty, deta OR a pierna 
en, E die? a A 


m 
SS 
3 

! 


ty”, + 2t cos(2r7ma)y; (2.47) 


m-—1 


EW = tm + 2tcos(27(m + 1)0)4m+1 (2.48) 


De acuerdo a las ecuaciones (2.47) y (2.48) a = p/q es un número racional, que 


permite calcular el espectro de energía del grafeno usando una red periódica de periodo 
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2q, siendo por tanto el problema de eigenvalores (2.47) y (2.48) expresado por una 
matriz de dimensión 2q x 2q dada en (2.49) con condiciones de frontera periódicas 
Voq = e y290p,. 


0 0 e 0 mn 2t cos(Ora) 
t 0 2tcos(21a) os 
0 2t cos(27ra) 0 e 0 An 0 
0 0 0 sa 0 2t cos(2qTa) 0 
0 0 0 2t cos(2qTa) 0 t 
2t cos(Ora) 0 0 de 0 0 
(2.49) 


En el límite de campos magnéticos débiles, la figura muestra claramente los 
niveles de Landau. Cuando el campo magnético es alto estos niveles de Landau 


presentan sub bandas [89]. 


Energía / t 





Figura 16: Mariposa de Hofstadter para el grafeno 
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Capítulo 3 


Transporte electrónico 


Nuestro propósito en el presente capítulo es describir las herramientas necesarias 
para describir el transporte electrónico en sistemas de baja dimensionalidad, conocidos 
como sistemas mesoscópicos, dichos sistemas son el puente entre los sistemas 
macroscópicos y la escala atómica. La escala mesoscópica tiene longitudes típicas que 
van desde decenas de nm hasta algunos yum, y es aquella en la que los electrones de un 
semiconductor o un metal ponen de manifiesto su comportamiento cuántico. En estos 
sistemas mesoscópicos el trasporte electrónico se encuentra influenciada por el tamaño 
del sistema L£, la longitud de onda de Fermi Ap, el camino libre medio elástico l., y la 


longitud de coherencia de fase L¿. 


La longitud de onda de Fermi Ap se encuentra relacionada con la energía 
cinética del electrón y por tanto depende de la temperatura del sistema. Cuando 
el sistema tiene dimensiones comparables a Ap la descripción clásica fracasa, siendo 
por tanto el comportamiento cuántico de las partículas muy importante. En sistemas 


semiconductoras Ap tiene ordenes de magnitud mayores que en un metal. 


El camino libre medio elástico es la distancia promedio que recorre un eléctron 
antes de ser dispersado elásticamente por un ión del sistema l, = vp7. (donde up es 
la velocidad de Fermi y 7, el tiempo medio entre colisiones elásticas). En metales l, 
es del orden de los 10nm. En cambio a temperaturas del orden de los mili-Kelvin en 


materiales semiconductores de alta pureza l. puede llegar hasta los 100/um. 


La longitud de coherencia de fase (L¿ ) es la distancia dentro de la cual la función de 


onda electrónica preserva su fase. Esta longitud es de vital importancia para definir el 


26 


transporte en sistemas mesoscópicos. Para un sistema de longitud finita L, el transporte 
es coherente siempre y cuando la longitud L sea menor que la longitud de coherencia 
de fase L¿ (L < Lg). 


Dentro del transporte coherente pueden distinguirse a la vez tres regímenes de 


transporte: 


Régimen Balístico: se da cuando el camino libre medio /, entre colisiones elásticas 
es mayor que L, Ap < L< le. 


Régimen Difusivo: en este régimen, el camino libre medio /, entre colisiones elásticas 
es menor que L, Ap < l,¿ < L. 


Régimen Localizado: se da cuando el camino libre medio /. entre colisiones elásticas 


es mucho menor que £, l, < Ap < L. 


Asímismo cuando Ap < l. < L¿ < L tenemos el transporte en régimen 


macroscópico. 


3.1. Formalismo de Landauer 


El formalismo de Landauer, se utiliza para describir el transporte eléctrico de 
nanoestructuras en el régimen balístico utilizando la conductancia eléctrica 23/40/41), 
expresado como: 

2e? 


G= FT(Br). (3.1) 


Siendo T (Ep) los coeficientes de transmisión cuántica correspondiente a todos los 


canales cuánticos permitidos en el nivel de energía de Fermi Ep. 


El sistema utilizado para el presente estudio se representa en la figura y se 
encuentra constituido por tres partes, el contacto izquierdo (L), el punto cuántico (D) 
y el contacto derecho (RR), los contactos conectan el punto cuántico a reservorios de 
electrones termalizados que ocupan los estados de los contactos. La conexión entre el 
contacto y los reservorios se consideran ideales (conductores perfectos), de modo que 
los electrones que vienen de los contactos y pasan a través del punto cuántico son 


dispersados solamente en la interfase entre los contactos y el punto cuántico. 
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nano 
dispositivo 


(D) 





Figura 17: Modelo de nanodispositivo conectado a dos contactos semi-infinitos 


La dimensión de los contactos es semi-infinita, de modo que dichos contactos 
son invariantes en la dirección de propagación y tienen infinitos estados permitidos. 
La dimensión del punto cuántico en la dirección transversal cuantiza los modos de 


propagación (P, permitidos, modos que son conocidos también como canales permitidos. 


Para la deducción de (3.1), se considera que la población de estados electrónicos 
en los contactos es dado por la distribución de Fermi-Dirac f(E — pa), donde pu, es el 
potencial químico del contacto a, representando a = L al contacto izquierdo y a = R 
al contacto derecho. La corriente eléctrica generada por los electrones del contacto Q 


que pasan a través del punto cuántico utilizando el canal p,, es dado por: 
2e 
Lan = 55 2 TEJA) FE=hHa), a=L,R. (3.2) 


donde W es la longitud del punto cuántico, el factor 2 es debido a la degeneración 
que presentan los niveles electrónicos por la presencia del espín electrónico, v(k) es 
la velocidad de grupo de los electrones y T,, (4) la probabilidad de transmisión de 


electrones por el canal (p, a la energía E(k). 


La corriente neta que pasa por cada canal p,, a través del punto cuántico se obtiene 
calculando la diferencia entre la corriente /7, producido por los electrones del contacto 


izquierdo y la corriente /r, producidos por los electrones del contacto derecho 
2e 
In = Min — Irnl = W Y u(k)T, (E) [F(E — uz) — FE — un)l. (3.3) 
k 


de modo que la corriente total a través del punto cuántico es la suma las corrientes /., 


sobre todos los canales permitidos: 


1= YI, = 2 DOTE) [HE — 1) — HE — un) (3.4) 


n,k 


en la ecuación ultima la sumatoria discreta se transforma al dominio continuo de k, 


2 > > Jar (3.5) 
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considerando la regla usual 


que permite reescribir (3.4) como 


1 ES farol TEE — mu) — HE — 1). (3.5) 
y considerando que 
_1dE 
v(k) = A dk” 
encontramos de 
IPN Jere HE = 14) HE ua). (37) 


a bajas temperaturas y pequeños voltajes aplicados eV = uz, — pg, podemos hacer una 
expansión en serie de Taylor de f(E — uu) — F(E — ua) al rededor de la energía de 
Fermi Ep, que conduce a f(E — u,) — F(E — ua) 2 Ó(E — Ep)(u — pa), resultado 
que permite expresar la Ec. como 


IPN JETS E — Ex) (ue — e). (38) 


considerando que la conductancia eléctrica G se define como G = I/V tenemos 


a= Y T,(8r) (3.9) 
=> h > n F ) " 
de modo que la fórmula de Landauer queda como: 
2e? 
G(Ep) = -T(Er) (3.10) 


siendo T(Ep) = > T, (Ep) el coeficiente de transmisión cuántica total sobre todos los 


canales permitidos. 


3.2. Formalismo de la función de Green 


Las funciones de Green son una herramienta poderosa para resolver muchos 
problemas de la física, en la presente sección describimos el método de la función de 
Green para calcular los coeficientes de transmisión, útiles para calcular las propiedades 


de transporte de sistemas nanoscópicos como el grafeno. 
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3.2.1. Función de Green 


En mecánica cuántica para describir una partícula sujeta a un potencial 


independiente del tiempo se resuelve la ecuación de Schódinger 


HAl|Yn) => Elm), (3.11) 


donde H es el operador hamiltoniano asociado a la partícula. Para resolver la ecuación 
(3.11) en representación (|r)) definimos la función de Green correspondiente como 
solución de la ecuación 


[E— H]G(r,r', E) =8(r—r'). (3.12) 


Si utilizamos el operador G(4) para representar la función de Green, entonces la Ec. 
(3.12) se reescribe como, 
[E — H]G(E) = 1, (3.13) 


de modo que 


G(r,r”, E) = (rIG(E)|r”. (3.14) 


Si E no es un eigenvalor de H, entonces la solución formal de la ecuación (3.13) para 


G(E) es dado por 
1 


E-H 
luego si E,, son los eigenvalores de H asociados a los autoestados normalizados |1,,) 


entonces se expresa como 
(9) yn 
G(E + 2 [bn) (| = 2 (3.16) 





G(E) = = [E -— H]”? (3:45) 





donde los E, son reales debido a que H es hermitiano. Por tanto G(L) es analítica en 
el plano complejo, con excepción a aquellos puntos en las cuales E toma valores reales 
E = E,,, puntos en las cuales G(4) presenta una discontinuidad (es decir el límite de 


G(E) cuando E se aproxima a E,, tiende al infinito). 


Sin e si no se especifica las condiciones de frontera, la inversa del operador 
E-—H en la Ec. no está definido, en dichos casos, con la finalidad de levantar la 
discontinuidad de e a ) asumimos que E = Ein con y >0, de modo que definimos 
las funciones de Green retardadas G*(E) y las funciones de Green avanzadas G” (E) 


como: 





G+(E)= lím G(E +in) (3.17) 


n>0+ 
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La Ec. (3.13) permite expresar G*(E) como 





G+(E) = a [E+in-H)* (3.18) 


en adelante nos referimos como funciones de Green G a las funciones de Green 
retardadas G*, de modo que G7 = (G+)f = GT. 


3.2.2. Funciones de Green fuera de equilibrio 


Las funciones de Green fuera de equilibrio NEGF son regularmente utilizados para 
calcular intensidad de corriente y densidad de carga en dispositivos a escala nanometrica 
bajo un cierto voltage aplicado. Este metodo es usado principalmente para el transporte 
balistico, pero puede ser extendido para incluir dispersion inelastica [42]. En esta 
sección derivamos la ecuación de NEGF para la corriente eléctrica y su relación con el 


coeficiente de transmisión. 


Siendo nuestro propósito determinar la función de Green Gp del nano-dispositivo 
de la figura [17| reescribimos la ecuación de Schródinger (3.11) para el sistema completo 


H, Hip 0 lL) [WL) 
Hp, Hp Hpr Ip) | =E | [Wp) (3.19) 
O Harp Hr [Yr) Vr) 


donde Hp es el hamiltoniano del nanodispositivo, H, y Hz los hamiltonianos de los 
contactos izquierdo y derecho, H,p con a = L,R son las matrices de acoplamiento 
entre el contacto a: = L, R y el nano-dispositivo, |Wp), WL) y (Wa) los estados cuánticos 


asociados a Hp, H, y Hp respectivamente. 


¿Por qué calculamos las funciones de Green? 


Para un sistema representado por el hamiltoniano A, las funciones de Green 
proporcionan la respuesta del sistema a una perturbación constante |v), en la ecuación 
de Schródinger: 

Hp) = Elp) + lv) (3.20) 


(E— H)l) = —[v) (3.21) 
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de modo que la respuesta a la perturbación es 
lp) = —[E — H]” |v) 


lv) =—G(E)|v) (3.22) 


La función de Green nos permite determinar muchas propiedades de un sistema 
físico, tales como la función espectral, la densidad de estados (DOS), la función de 
transmisión, etc. Por ejemplo calculamos los estados |, ) asociado al contacto izquierdo 
y l4) asociado al contacto derecho en función de los estados del punto cuántico |4p) 
a partir de la ecuación de Schródinger (3.19), que se reescribe como: 


Hrlvi) + Hipi4p) = Elvr) (3.23) 
Hpilvz) + Holvp) + HorlVa) = Elvp) (3.24) 
Hraolpo) + Halva) = Elvr) (3.25) 


obteniendo |4,) de la Ec. (3.23) y lr) de (3.25) tenemos 
1) = [E — Hi] Hiolbo) = 91H LD |Do), (3.26) 


lr) = [E — Ha] Hrolp) = grHrolWo). (3.27) 


donde yg, =|E — H¿]* y gr = [E — Ha]? son las funciones de Green asociados a los 


contactos aislados izquierdo y derecho respectivamente. 


Auto-energías 


Para obtener la función de Green retardada Gp asociada al nano-dispositivo de la 
figura [17] hacemos e = E+1n y reescribimos la definición de la función de Green (3.13) 


correspondiente a la ecuación de Schródinger es (3.19) como 


e=- Hi —H¡p 0 GL Gp Gr IT00 
—Hpz e= Hp —Hpr Gp Gp GprR = 010 (3.28) 
0 —Hrp e=- Hr GrL Grp Gr 0.017 
seleccionando tres ecuaciones de (3.28) tenemos 
(e — H¡)Gip — HipGp =0 (3.29) 
—Hp1G in + (e- Hp)Gp — HprGrp =] (3.30) 
—HrpGp + le = HRGrD =0 (3.91) 
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obteniendo Gp de y Grp de y reemplazando en (3.30) se tiene 
[—Hpile — Hi) “Hip + (e - Hp) —- Hprle- Hi) *Hap] Gp =1 (3:32) 
considerando las auto-energías asociadas a lo contactos como 
2, = Hpigi Hip (3.33) 


2 = Hpr9rHArD (3.34) 


donde g, = (e— Hi)? y gr = (e— Hp)? son las funciones de Green correspondientes 


a los contactos izquierdo y derecho aislados, entonces la ecuación (3.32) se escribe como 
[—Es + [e— Hp) - 3aJGp =1, (3.39) 

luego la función de Green del nano-dispositivo Gp es 
Gp(E) = (el - Hp -YE,-—YEn)”* (3.36) 


Las matrices de acoplamiento ',(E) y Pr(E) se definen como la diferencia entre las 


auto-energías avanzadas y retardadas de los correspondientes contactos 
Tor = (E, — El p) (3.37) 


De acuerdo a la ecuación (13.36) podemos decir que el efecto de los contactos en el 
dispositivo es adicionar las autoenergías al Hamiltoniano del dispositivo, es decir las 
autoenergías transforman el sistema completo de la figura [l8fa en un nuevo sistema 
formado por el dispositivo mas autoenergías representado en la figura[18fb. En términos 
de cálculo computacional en lugar de calcular la función Green del hamiltoniano 
asociado al sistema completo, calcularemos la funcion de Green para el hamiltoniano 
efectivo H.¿ = Hp +*1 + *n. Sin embargo debemos tener en cuenta que este es un 
proceso para calcular la funcion de Green, dado que los autovalores y los autovectores 


de este Hamiltoniano efectivo no son cantidades que se puedan interpretar facilmente. 
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contacto contacto 
izquierdo i iti derecho 















nano dispositivo + autoenergías 


(b) 


Figura 18: Las autoenergías permite reducir el sistema completo (a) de dimensión infinita 


a un sistema finito (a) de dimension igual a la del nano-dispositivo. 


3.2.3. Coeficiente de transmisión y la función de Green 


Deducimos la expresión que nos permite calcular el coeficiente de transmisión en 
función de las funciones de Green, para ello evaluamos primero la respuesta que produce 
un nano-dispositivo a las ondas incidentes, considerando que los electrones (ondas 
incidentes) que se inyectan en el nano-dispositivo ocupan los estados correspondientes 


a los contactos [44]. 


Sea |4z n) un estado que corresponde al contacto L que para cierta energía E es 
reflejado completamente en la interface entre el contacto y el punto cuántico con n un 
modo de propagación y |y+) el estado retardado obtenido como respuesta del sistema 


completo, entonces el estado total | n) + |1W*) satisface la ecuación de Schródinger 
H(ben) + 107) = Elld1m) + 117) (3.38) 


Albin) + Hlw*) = Ellúon) + l1*)) (3.39) 
(Hp + A+ Hipo + Hp, + Hor + Hao + Ar)ldon) = (AL + Hop) (rn) 
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debido a que |Wz .) es una estado que corresponde solo al contacto L 


luego 


(A, + Hor)lWin) + H iy?) = Elli) + 1*)) 
(E+HoL)lWin) + Hly*) = Elvin) + Ely?) 
Horlbin) + Hly*) = Ely?) 

(E pe H) y?) == HbLlWLn) 








ly+) =[E— H]*EpilWin) = GHorlzn) (3.40) 


Con este resultado el estado del punto cuántico |4p) se escribe como 


lp) = GoHorlVn). (3.41) 


Además los estados de los contactos (3.26) y (3.27) se reescriben como: 


lv) = (1+9,Ao0GpHpr) Win) (3.42) 
[Ve) IrRArDGDA pi lULn) (3.43) 


La corriente eléctrica inyectada por el contacto L y pasa al contacto R a través del 


punto cuántico es dado en función de |4p) y |Wp) como: 


Laca = 5 nl Hol) — (vol Hontn) (3.44) 


reemplazando |4p) y |) dadas por (83.41) y (3.43) se tiene 


Lo41 


le 

=- ((WalHoG por.) = (nl Ho Gh Honlbr)) 
le 

==. (We oCh Hor HroGoHorlbrn) ES (Won HioGhHorgrHroGo Hor Ven) 
le 

E (Wen Ho GLERGo Horn) — (Gen Mio GLExGD Horn) 


(Won HipGhi(Er y DG pH) 


>] o 


e 
Las 5 rl Hip CRC DH pls). (3.45) 


La corriente total que pasa del contacto izquierdo al derecho a través de todos los 
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modos n, considerando la degeneración de espín es 











ba = 7 a dEF(E — pu) Y AE — En HioGhTrGoH 11) 
= Ñ dEF(E — ur) Y (E — En Mora] Hip lm)(m¡G TG pH pr) 
= 7 Y dEf(E — uy) y (E — En (miG5TrGoHoilWin (don Hi p|m) 
= a Ñ dEF(E—yy) za GhTrGpHpr (y ó(E — Eeualtóa) Hp|m) 
== a Ñ dEf(E — yy) 20 Gh rGpHpr (9 E 3.) H pim) 
= o Ñ dEf(E—p) 20m GpIrGp (Ez e E) |[m) 
= o Ñ dEf(E—p) 2 Gp rGp! L|m) 
la = = 1 _ dEF(E — y Tr[G TaGpT7). (3.46) 


De manera equivalente se encuentra que la corriente total que pasa del contacto derecho 
al contacto izquierdo es 
O 


hiea==>>2 |] dENÑE- paTriG TI rGpT1). (3.47) 
La corriente total a través del punto cuántico es entonces: 
le f% 
Ia) CEE) HEM) TAG AGO] (848) 


comparando las corrientes en el punto cuántico dadas por las Ecs. (3.7) y (3.48) se 
tiene que 

T(E) = Tr[G TG pT1). (3.49) 
La Ec. (13.49) es equivalente a la relación de Fisher-Lee que expresa el coeficiente 


de transmisión en términos de las matrices de transmisión como: 


F=THEL (3.50) 


3.2.4. Método de decimación 


Utilizado para determinar las funciones 9, 9r y Gp de manera recursiva [47], debido 


a que los contactos se consideran de dimensión semi-infinita, y el punto cuántico tiene 
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dimensión finita. La aplicabilidad del método de recurrencia se debe a que tanto los 
contactos y el nano-dispositivo se pueden descomponer en sub redes verticales finitas 
en serie como muestra la figura donde cada sub red es llamada en adelante como 
cadena unidad, dentro del modelo del tight-binding cada cadena unidad interactúa 


solamente con las cadenas unidad laterales próximas. 


“=. -4 -3-2-1 0.1 2 3 4... 





Figura 19: Aproximacion de un nanodispositivo utilizando la cadena unidad. 


Para la presente deducción consideramos un sistema compuesto por cinco cadenas 
unidad, que en la figura se denotan por —2,—1,0,1,2, siendo el hamiltoniano 
tight-binding asociado dado por 


Ho Hoi 0 0.0 


Ho, Hai Hao 0 0 
H= 0 Hio Ho Hor 0 (3.51) 
0 0 Hi, Hi Hi, 
0 0 O Hi, Ha 
donde H,,, con n = —-2,-—1,0,1,2 es el hamiltoniano correspondiente a la n-ésima 


cadena unidad, H,, +1 es la matriz de acoplamiento de la cadena n con la cadena n+1 
y Elsa es la matriz de acoplamiento de la cadena n + 1 con la cadena n. La función 
de Green asociado al hamiltoniano (3.51) se define de acuerdo a (3.13) por 


a o o a Ga Gai Gao Gar Cas 1.0.0 0.0 
Ha € H-1  —H-io0 0 0 Gia  G_i Gio CGií1 Gia 0. 1.0.0.0 
0 Hijo <-Ho —Ho1 0 Go-2  Go-1 Go Go1  Goz |=|o o 1 0 0 
0 0 Hb e= Hi —-H1 2 G1,-2 G1,=1 Gi,0 Gi G1,2 0.00. 0 1.0 
0 0 0 Hi, «e-H3) NG2-2  G2a-1  Gzo  Ga21 Ga 0.0.0.0 1 


(3.52) 
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esta ecuación matricial se reescribe como 


le=B9)6 ap = Ba aGas=0 (3.03) 

—H, ¡Goo + (e- H_1)G_10— H_10G09=0 (3.54) 
—H 1 0G-10 + (e— Ho)Go — Ho1Gio =1 (3.55) 
—HÍ,Go + (e — Hi)Gio0 — Hi 2Goo =0 (3.56) 

—H| ¿Gio + (e — H3)Go,o =0 (3.57) 


De la Ec. obtenemos la función G20 
Gao = E H>H ¿Gro 
y sustituimos en la Ec. 
—HÍ,Go + (e— Hi)Gio — Hio(e — Ha) *HÍ ¿Gio =0 
obteniendo G¡ y tenemos 
Gio =(e— H¡ — Hio(e— Ha) *H] 2) *H],Go (3.58) 

Asimismo de la ecuación (3.53) obtenemos la función G_20 

Gu =18=H 2) *H 39-16-10 


y sustituimos en la Ec. (3.54) 


—H,_,(6 — Ha) *H_3-1G-10+(e- H_1)G_10- H_10Go =0 








obteniendo G_¡ y tenemos 
G_10 = le = H_; = Hoale — H y E 2, y Ha 1LoGo (3.59) 
sustituyendo (3.58) y (3.59) en la Ec. (3.55) tenemos 


Hole Hg HU le — H_2)*H_2-1)*H_10Go + (e — Ho)Go 
—Horle — Hi — Hio(e — Ha) *Hl))*H] ¡Go =0 








de aquí obtenemos que la función de Green Gy se expresa como 


=1 








Go= |e-H,-H,) Hi Hoy Hi 
e— H, — ES — 








(3.60) 
Para el cálculo de Gp esta recursión se extiende hasta el tamaño del punto cuántico, y 


en el cálculo de g, y gn se considera hasta un tamaño suficientemente grande. 
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3.2.5. Funciones de Green de contacto derecho 


En la bibliografía existen varios métodos recursivos para determinar las funciones 
de Green de los contactos, así tenemos el método auto-consistente [46], el método de 
decimación [23/lu8l[49]. Nosotros utilizamos una variante del método de la matriz de 
transferencia [501152], para calcular de manera rápida y eficiente la función de Green 
del contacto derecho, la particularidad de este método es que se aplica cuando los 
hamiltonianos asociados a cada una de las cadenas unidad del contacto derecho son 
todas iguale. Para el presente cálculo consideramos que el contacto derecho es una red 
semi-infinita, que de acuerdo a la figura es representado por las cadenas unidad 
n = 0,1,2,... de modo que la ecuación que define la función de Green del 


contacto derecho se escribe como: 


e=- Ho  —Ho1 0 0 pe Go Go Go Gos 100.0 
-H, «-Hy -—Hp 0  ...||Gw Gu Gir Gi 01.00 
0 —Hl, e=- Ho —Hoa3  ... Go Goa Gro Ga =|001 0 

0 0 —Hl, E= H33 a G30 G31 G32 G33 0.00 1 
(3.61) 


considerando que Ho = Han+1, Hio = Hn+1m y Hnm = Hoo para todo n tenemos 





(€ — Ho0)Goo — HoGio = 1 (3.62) 

—Hi0Goo + (e — Hoo)Gro — HoGao = 0 (3.63) 
—Hi0Giro + [e — Ho0)G20 — HoGzo = 0 (3.64) 
—Hi0Gao + [e — Ho0)G3o — HorGao = 0 (3.65) 
—Hi0G3o + (e — Hoo)Gao — Ho1G5o = 0 (3.66) 
—Hi0Gao + [e — Ho0)G50 — Ho1G6o = 0 (3.67) 
(3.68) 


despejando Go de ((3.63) tenemos 


Gio = (e Hoo) (HioGoo + Ho1G»o) (3.69) 
reemplazando (3.69) en (3.62) encontramos Goo 
(€ — Hoo)Goo — Horle — Hoo) (HioGoo + Ho1Gao) = 1 
(e = Ho = Horte => Hoo) *H10)Goo = I+ Hor (e = Hoo) *Ho1Goo (3.70) 
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Despejando ahora Gzp de (3.65) se obtiene 
G3o = (e — Hoo) *(Hi0G20 + Ho1Gao) (3.11) 


reemplazando (3.69) y (3.71) en (3.64) tenemos 


—Hio[(e — Hoo)*(Hi0Goo + Ho1Ga0)] + [e — Hoo)Gao — Hor[(e — Hoo) *(Hio0G2o0 + Ho1Gao) = 0 


reordenando 


[e — Hoo — Hiole — Hoo) Ho — Hule — Hoo0)”*Hi0o]G20 


(3.72) 
= Hiole — Hoo) *Hi0Goo + Hnle — Hoo) *Hi0Gao 
haciendo el cambio de variables 
H, = Ho+Hnle— Hoo)*Hio 
Bi = Hiol(e— Ho) *Hio 
o = Houle- Ho) *Ho 
H, = Hg + Hole — Hoo) *Ho1 + Horle — Hoo) Ho 
las ecuaciones (3.70) y (3.72) se reescriben como 
(e- Hi)Go = I1+01Go0 (3.73) 
(e Hi)Gzo = P1Go +01Gs0 (3.74) 


si despejamos Gso de (3.67) y tomamos Go de (3.71) para reemplazar en (3.66) y seguir 
el mismo procedimiento seguido para deducir (3.74) encontramos 


(e- Hi)Gao = P1G2+01G60 (3.75) 
así mismo encontramos que 


3.76 
3.17 
3.18 
3.79 


(le - Hi)Góo = P1Ga10+01Gso 
(e- Hi)Gso = P1G60+01G10 
(e — H¡)Gi0 P1Gs0 + 01612 


(3.76) 
(3.77) 
(3.78) 
(3.79) 


Para hacer una segunda iteración despejamos Gap de (3.74) y reemplazamos en (3.73) 


obteniendo 
(e- Hi)Go = I+o1((e— Hi) *(8B1Goo + 041Ga40)) 
le == H; = ale e Hi *B1)Goo = I + ale 7 Hi) *01Gao (3.80) 
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despejando Gap de ((3.74) y Góo de (3.76) y reemplazando ese obtiene 
(e - Hi)Gao = Brlle— Hi) *(B1Go0 + 041G40)) + a1((e — Hi) *(81Gao + 01Gsp)) 


(e- Hi — Bile—- Hi) 01 — ar(e— Hi) *B1)Gao = Bite — Hi) *81Goo0 + or (e — Hi) *01G48.81) 
haciendo el cambio de variable 
Hz Hi + ou(e- Hi) "81 
Ba = Bile-H) "Br 
om = o(e-H¡) aj 
EL = Hi+ ble HE) + ole HYTA 
escribimos y como 
(e—- H2)G00 = I+07Gg (3.82) 
(e — H3)JGgo = B2Go + 09Gg0 (3.83) 


de manera similar encontramos 


(e- Ho)Gso = P2Gao + 09G12 
(e — H,)G 12 = P2Ga40 + 09612 
(e- H3)Gió6 = P2Ggo + 0916 


De acuerdo a las ecuaciones de la primera iteración al y las ecuaciones de 
la segunda iteración (3.80) al (3.86), después de la n-ésima iteración tendremos 

3.88 
3.89 
3.90 
3.91 
3.92 


(e = Hn)Go = 1+0nGono 
(e —= Hi )Gono = BnGoo + AnGaran o 
(e — Hi, )Goxano = B2Goano + 02G3s2m 0 
[E Hr) Gaxon.o =  P2Gison o + 09 Garan o 
es Hr) Gaxon o =  B2aGaxon o + 095x920 





( 
( 
( 
( 
( 
(3.93 


) 
) 
) 
) 
) 
) 


con 


Ha — Hai ER Qn—1[€ + 1 
Bn => Bn-1[€ _ UA e sE 
dr == ales a a 


n 


Ha = a ole Hn-1) “On1 + Onil[e= Ha 


3.94 
3.95 
3.96 
3.97 


COS CTS CTS OS 
SU E a E E 
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Si el número de iteraciones n es suficientemente grande entonces se encuentra que 
0%, es considerablemente pequeña, de modo que la función de Green aproximada del 


contacto derecho es dado por: 


Goo == (e — Hn) 7. (3.98) 


3.3. Formalismo de Buttiker para 4 contactos 


Buttiker realiza una extensión al formalismo de Landauer que describe el transporte 
electrónico de nano-dispositivos con dos contactos a sistemas cuánticos con cuatro 
terminales [531154], en general este tratamiento se puede extender a más terminales. En 
principio consideramos un sistema multiterminal, donde cada contacto a tiene varios 
canales permitidos M,,, entonces la ecuación (3.7) para corriente total [, en el contacto 
o, se generaliza como: 

26 


Ta = dE 


na (Mo — Roa) fa — > Tap 


Bra 


(3.99) 








el factor 2 se debe a la degeneración del espín, si expandimos fg y f. al rededor del 


nivel de Fermi como 


falo) = 0) + (3g) Va 


considerando la aproximación hasta el primer orden, la ecuación (13.99) se escribe como 





La = 7 dE (Ma — Rao) (Luo) + (56 Eje va) - Ea (10) )+ (55) 1) 
a = 7 dE |fíup) (o e - Eno) + (Ma — Roa) (55) ev. a 





La conservación de corriente en el contacto a y la simetría de la conducción eléctrica 


nos permite escribir: 


(Ma — Rua) — Y Tap =0, (Ma — Roa) — Y Tgo =0 (3.100) 
pFa BHO 
luego /,, es 
le df df 
L,=3 | dE |(M. — Raa) (35) eV, — 2 Tos (ds) eVa (3.101) 
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asumiendo que 
Je? 2e? 
Gua sé Sr dE (e) (Ma E > Gas E E dE (ds) Tag (3.102) 


tenemos 


a N GagVa. (3.103) 


a bajas temperaturas Gua y Gag dados por (3.102) se expresa como 


2e? 2e? 
Gua = Ma NN Raa); Gas = Tap (3.104) 


Para un punto cuántico con cuatro contactos a bajas temperaturas la corriente /., 
en el contacto a (a. = 1,2,3,4) dado por (3.103) es 
2e? 
La = Y |(Ma — Raa)Va — Y TagVa|, (3.105) 
h 
pFa 
cabe recalcar en esta parte que M,, es el número de canales permitidos a través del 
contacto q, es decir que Roa = 7 R?, es la probabilidad total de reflexión de las ondas 
incidentes por los M,, canales permitidos en el contacto a, así mismo Tag = > Tag es 
el coeficiente de transmisión total de las ondas que inciden por el contacto q y pasan 
al contacto 6. 





Figura 20: Punto cuántico con cuatro contactos. 


Considerando que toda la corriente que ingresa a través del contacto 1 pasa al 


contacto 3 y toda la corriente que ingresa por el contacto 2 pasa al 4, es decir If, = —l3 
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y Ll, = —Í, a partir de (3.105) encontramos que 


Tia + T32 Ta + Ta 
S S 


S = Tia + Tia + Tao + Tia = Ta + Tar + Tos + Tas 


61. (14 =Yj= 1M=18) 


con estas identidades reescribimos la ecuación (3.105) como 


Li _2e? Yi. —Y12 V, — V3 
L> h A-v21 Ya Va — Va 











donde 
T TuNT T 
vi = (Mi, — Ri) - Ca + al 12 + Tia) 
e Ti2T34 — T3oTia 
12 S 
_- TaTag—TaTa 
ya = S 


(Ta + Ta3) (To + Tia) 


Ya — (Ma ES Raz) > S 





(3.106) 


(3.107) 


(3.108) 


(3.109) 
(3.110) 
(3.111) 


(3.112) 


para la medida de la resistencia Hall, consideramos /2 = 0 en los contactos 2 y 4, de 


modo que (3.108) se reescribe como 


des 

L = q aulV NN V3) => al Va = Va)] 
de? 

0 = y Hra(Va — V3) + yaa (Va — Va)] 


resolviendo para /, encontramos 











Ze? De? = 
ha == a (de Va) We Va) ES 11722 A Va) 
h 21 h 21 
entonces la resistencia Hall Ry es 
Va — Va h y21 
Ra =R = = 3.113 
di dd Í¡ 2e? Y11Y22 — Y21Y12 ) 
y la resistencia longitudinal Ry es 
Vi, — Va h 
Aya HA 122 (3.114) 
Li Zes Yes — Ya Ya 


En las figuras y [23] presentamos los resultados de Ry, Ry y DOS para un 


punto cuántico de grafeno formado por una red finita de dimensiones 23a x Ly 3a, 
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siendo el ancho de los contactos izquierdo (1) y derecho (3) igual a 23a y el ancho de 
los contactos superior (2) e inferior (4) igual a 4 Y3a. El campo magnético uniforme 
aplicado sobre el sistema es representado por el potencial vector A = — Byy2 y su efecto 
sobre el punto cuántico es dado por la sustitución de Pierls. Para fines de obtener los 
resultados calculamos primero las funciones de Green de los contactos por medio del 
método recursivo rápido descrito en luego calculamos las auto-energías y las 
acoplamos al punto cuántico siguiendo las matrices de transformación propuestas por 
Jayasekera [55], seguidamente determinamos la función de Green, finalmente utilizando 
la función de Green calculamos la resistencia Hall Ry por medio de la ecuación (3.113), 
la resistencia longitudinal R, usando la ecuación y DOS respectivamente, el 
programa numérico utilizado para estos cálculos es la que proponemos en [56]. Otra 
alternativa de tratamiento de las auto-energías para sistema nanoscópicos mas grandes 
es utilizar las propuestas de Thorgilsson [57]. 


La figura representa el diagrama de colores de la resistencia Hall obtenida a 
partir de en la gráfica el eje horizontal representa el flujo magnético por unidad 
de flujo cuántico y el eje vertical la energía (E /t), los colores cálidos representan los 
valores positivos que toma Ry y los colores frios los valores negativos de Ry, la gráfica 
muestra también la naturaleza anti-simétrica de Ry respecto a E = 0. Así mismo 
observamos la formación de los niveles de Landau, niveles que se aprecian mejor en la 
gráfica de la DOS dada en la figura [23| y en la resistencia longitudinal R, dado por 
(3.114) y representado en la figura [22] La presencia de los niveles de Landau se deben 
a la localización de los estados de borde en determinados valores de energía a medida 
que el aumenta el campo magnético y son indicadores que prueban la cuantización de 
Ry en el punto cuántico estudiado. La cuantización de Ry es mas notoria a medida 
que las dimensiones del punto cuántico se incrementan, ya que cuando las dimensiones 
del punto cuántico son muy pequeñas los efectos de cuantización espacial son mas 


preponderantes. 
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0 0.1 0.2 0.3 0.4 
a = pig 


Figura 21: Resistencia Hall Ry dado por la ecuación (3.113) para un punto cuántico de 


Grafeno con cuatro contactos 


En la figura [22] los colores cálidos representan valores positivos de Ry y los colores 
fríos los valores negativos de Ry, así mismo se observa que Ry tiene propiedades 
simétricas respecto a E = 0 a diferencia de Ry. La figura [23] es el diagrama de colores 
de la densidad de estados del punto cuántico que muestra claramente la presencia de los 


niveles de Landau y la concentración de los estados de borde en los niveles de Landau. 
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Figura 22: Resistencia longitudinal Rz dado por la ecuación (3.114) para un punto 


cuántico de Grafeno con cuatro contactos 
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Figura 23: Densidad de estados de punto cuántico de Grafeno con cuatro contactos 
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Capítulo 4 


Nanoribbons de grafeno 


Los nanoribbons de grafeno son bandas (cintas) de grafeno con anchos muy 
pequeños (< 50 nm) como muestra la figura En estos sistemas la distribución de 
los átomos en los bordes del nanoribbons juegan un rol importante en la caraterización 
de sus propiedades electrónicas. En el presente estudio clasificamos lor bordes de 
los nanoribbons de grafeno en bordes zigzag, armchair y bearded (barbado) como se 
muestra en la figura [19][27][58][59]. A continuación describimos analíticamente la 
estructura electrónica y los estados electrónicos de nanoribbons de borde zigzag y 
armchair utilizando la ecuación de Dirac y el modelo de tight-binding, calculamos 
la conductancia eléctrica y analizamos el efecto del campo magnético sobre los 


nanoribbons de diferentes bordes. 


1) 
. 
1) 
.o .o e... pe. 





.s .o o... 


Figura 24: Nanoribbons de borde zigzag con ancho finito en el eje y. 
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Figura 25: Red finita de grafeno que muestra los tipos de borde zigzag, armchair y 


bearded utilizados en la descripción de los nanoribbons. 


Para obtener la estructura de bandas y eigen-funciones de los nanoribbons 
reescribimos el hamiltoniano (2.33) en términos del vector de onda k = (kz, ky) al 
rededor de los puntos K y K” como 


0 ko—ik, 0 0 

H= hop | He +0 0 o di 
0 0 0 — (Ko + ikiy) 
0 O —(ko—ik,) 0 


Las soluciones a la ecuación de Dirac para el hamiltoniano (4.1) se expresan como 


my = [Ux(1) 
V(r) (0) (4.2) 


donde Yx y Yx, son los espinores asociados a los puntos K y K”, las funciones Vx y Uk: 
tienen cada una dos componentes denotados por Y%, vE y Y, 1112 respectivamente. 
Los superindices A y B en las componentes de los espinores representan a las subredes 


A y B del grafeno, de modo que la función de onda total de la nanocinta se expresa 


os p*(r) — ¿Kr lr) ¿Kr Ve (r) 
n= (0) (io) * o) je 


en el tratamiento siguiente consideramos que el ancho de la nanocinta se encuentra en 


como 


el eje y y su longitud considerada infinita se encuentra a lo largo del eje 1, de modo 
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que debido a la invariancia traslacional en el eje x las funciones de onda deben cumplir 


la condición de Bloch 


v(r) = (y). (4.4) 


4.1. Estructura electrónica de nanoribbons de grafeno 


Z18Zag 


Analizamos la estructura electrónica del nanoribbons de borde zigzag de ancho 
W = (2 — ) a, con N número de dímeros en la cadena unidad mostrado en la figura 
Haciendo que ik, —> 0, la ecuación de Dirac para el hamiltoniano (4.1) se escribe 


como 


0 kz — Oy 0 0 vi ve 
k B B 
e sa p i a 0 (4.5) 
0 0 0 — (Ko + 0.) Vir Vir 
0 —(k => 0,) 0 VE VE 


de (4.5) las ecuaciones para evaluar 44 y vé en el punto K son 


a A 


(ko — Oy) VR = EVE 


4.7 
(ha + 0,)ud = EUR ll 
con £€ = €/hvp, desacoplando las ecuaciones 
1 
Vx = ¿(ho — Oy) g (4.8) 
1 B_ 2B 
(ko Ap 0y) 3 (Ko = O) Vk = EVk 
(Eo — OE = E PK 
OjOz + (E — k)0z =0, (4.9) 
la solución general de la ecuación (4.9) es 
vrly) = ae + Pen (4.10) 
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con y = y//k2— £2 que puede ser real o imaginario. Reemplazando 4% (y) en (4.7), se 
obtiene que (y) es 


vély) = (ke — Ne” + Blka + ye" (4.11) 


Debido a que el ancho del nanoribbon en el eje y es finita, asumimos que los bordes 


están dados por y =0 y y = W, entonces las funciones de onda deben satisfacer 
$ (y=0)=0, dy =W) =0 (4.12) 


La condición ¿%(y = 0) = 0 se asume debido al hecho de la frontera y = 0 del 
nanoribbon está constituida por átomos de la subred B y la condición 9*(y = W) =0 
debido a que la frontera y = W del nanoribbon es constituida por átomos de la subred 


A. la aplicación de estas condiciones de frontera a la Ec. (4.3) conduce a la expresión 


o (A =00Y _ a (Ra (UR =0 
9% (y =W) Wiely = W) biely = W) 

0=e EE (y =0) + =p, (y =0) 

0=e Ely =W) + pg (y = W) 


condiciones que se verifican si 
vily =0) = vi ly =0) = vi ly =W) = vi ly =W) =0 (4.13) 
haciendo que verifique 42 (y = 0) =0 se tiene 
pE(y=0)=0=a0e0 4 Be10 


a=-8 
luego usando 4 (y = W) =0 en (4.11) tenemos 
vily =W) =0= a(ka — ye" + Bllo + ye 
0=a(k, —- ye” — 0(k, + yjeW 


quedando la ecuación trascendental 


a (4.14) 





La ecuación tiene soluciones reales para y si k, > 1/W, soluciones que 
corresponden a funciones de estado de borde (edge states) localizados en la frontera 
del nanoribbon, con energías 

E =HvY ki — y, (4.15) 


el resultado (4.15) no puede ser predicho a partir de la estructura de bandas del grafeno. 


En otro caso si k, < 1/W las soluciones de la Ec. (4.14) son puramente complejas 


de la forma y = ik,,, entonces 


kg — Mn —2ikn W 
e 1 = e 
Kg + ikn 
(ko —ikpJe e" = (kg + ikpJe tr W 
Eden Em mW) pe ik [enW de gn 


kg sin(kW) = kh, cos(k,W) 


quedando 


kr, 
O A 4.1 
Ko tan(k, W) o 


las soluciones de (4.16) corresponden a estados confinados en la nanocinta, que tienen 


E=+ YM + (4.17) 


energías 


4.1.1. Método del tight-binding 


La dispersión de energía de los nanoribbons dentro del modelo del tight-binding 
describe de manera eficiente su estructura de bandas. El hamiltoniano asociado a 
los nanoribbons se escribe de manera equivalente al hamiltoniano del grafeno, 
utilizando como condición de frontera que las funciones de onda en los bordes del 
nanoribbons se anulan, como en la ecuación (4.12). De acuerdo a la idea propuesta 
por Lee y Joannopoulos para el cálculo de funciones de Green de conductores 
semi-infinitos, consideramos que el nanoribbons se encuentra formada por cadenas 
verticales finitas, siendo la cadena unidad del nanoribbons zigzag representada dentro 
del rectangulo rojo en la figura cada cadena unidad se repite periódicamente a lo 
largo del eje del nanoribbons (eje 1) e interactúa con las cadenas laterales y vecinas a la 
izquierda y derecha, de modo que el hamiltoniano del nanoribbons N-zgnr es expresado 


por la matriz de dimensión 2NV x 2N siguiente 
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0 2t cos(k.a) 0 0 0 0 
2t cos(k.,a) 0 É 0 0 0 
0 t 0 2t cos(k,a) 0 0 
H= 0 0 2t cos(k,a) 0 t 0 
0 0 0 $ 0 2t cos(k.a) 
0 0 


0 0 2t cos(k.a) 0 
(4.18) 
La estructura de bandas del nanoribbons zigzag calculadas diagonalizando el 
hamiltoniano (4.18) se ilustra en la figura|26|para seis diferentes anchos del nanoribbons. 
Los resultados muestran que la cantidad de estados permitidos para cada energía 


depende directamente del ancho de los nanoribbons. Otro aspecto importante que 


se observa en los nanoribbons de borde zigzag es su comportamiento metálico 


independientemente del ancho del nanoribbons [351621163]. 
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Figura 26: Espectro de energía de nanoribbons de diferentes anchos (a) 7-zgnr, (b) 8-zgnr, 
(c) 9-zgnr, (d) 10-zgnr, (e) 30-zgnr y (f) 100-zgnr 


4.2. Estructura electrónica de nanoribbons Armchair 


Analizamos nanoribbons armchair de ancho W = (N-—1) a, donde N es 
el número de dímeros en la cadena unidad mostrado en la figura así mismo 
consideramos que los bordes del nanoribbons son dados por y = 0 y y = W, en 
estas fronteras de acuerdo a rotación en 7/2 de la figura [3| tendremos átomos de las 


subredes A y B simultáneamente. 
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Figura 27: Nanoribbons de borde armchair con ancho en la dirección del eje y. 


La deducción del hamiltoniano de Dirac (4.1) se hizo considerando que el borde del 
nanoribbon era la dirección zigzag, entonces para el estudio de nanoribbons armchair 
introducimos una rotación de 7/2 en el vector de onda, es decir hacemos k.,, —> —ky y 


ky > k, luego si consideramos que ky —> 10, entonces el hamiltoniano (4.1) se escribe 


como: 
0 iO, + ik 0 0 we, Dé 

hup có : UR =8 UR (4.19) 
0 0 10, dez 0 pe vE 


de (4.19) las ecuaciones para obtener 4% y YE en el punto K son 


—i(O, + ko) VR = EVK 


Ñ (4.20) 
(0 ka) = EVR 
desacoplando las ecuaciones 
(0) + ko = EV (4.21) 
con y = A, B. Las soluciones para 4% y vÉ son 
vA(y) = Actry 4 Bei, (4.22) 
hey) =CO "EF Den (4.23) 


como en los bordes del nanoribbon existen átomos del tipo A y B, las funciones p% y 
6% definidas por (4.3) deben cumplir las condiciones 


entonces tenemos 





0=rpily =0) + epa (y =0) (4.26) 
0= "pz (y =0) + pg, (y = 0) (4.27) 
p= AE W ¿mat (y = W) + e EW ¿mae (y = W) (4.28) 
0=¿W¿mtgEly =W) + e EW ¿mer (y =W) (4.29) 
estas condiciones se verifican si 
Vily =0) + Y ly =0) =0 (4.30) 
My = W) + EM (y =W) =0 (4.31) 
Reemplazando (4.22) y (4.23) en estas condiciones tenemos 
A+4B+C+D=0 (4.32) 
AQEimtEIW y pam EW y O JW y DW — (y (4.33) 
escogiendo B = C =0, se tiene A=-—D y 
A Ca pa En 
sin[(y, + K)W] = 
(yn + JW =nwT 
NT K= nT 4 (4.34) 
Vn W W  3v3a . 


Reemplazando y = AeY en la Ec. (4.21), se encuentra 
Oy Vx = iynAe?, Ox = (iy) Ac = — 
(rm) + ha) 0% = Evk 
mt = E 


¿=tyYN+q (4.35) 
dado que todas las soluciones de y, conducen a soluciones oscilatorias, se concluye que 


no existen estados de borde localizados. 


4.2.1. Método del tight-binding 


Describimos la dispersión de energía de nanoribbons de borde armchair usando el 
hamiltoniano tight-binding dado para el grafeno, utilizando de acuerdo a la figura 
[27] las condiciones de frontera dadas en la ecuación y que establecen que 
las funciones de onda de las subredes A y B en los bordes del nanoribbons se anulan 
simultáneamente. Así mismo consideramos que cada cadena unidad representada por el 
rectangulo rojo en la figura[27]se repite periódicamente a lo largo del eje del nanoribbons 
(eje 1), cada cadena unidad interactúa con las cadenas laterales y vecinas a la izquierda 
y derecha de modo que el hamiltoniano del nanoribbons N-agnr es representado por 


una matriz de dimensión 2NV x 2N expresada como: 


0 t 0 0 terre 0 0 
t 0 t 0 0. teta q 
0 t 0 t 0 O tela 
0 0 t 0.0 0 0 
H= | gia 0. 0.0 t 0 06) 
0 tera 0 GU 0 t 
0 DO tee 0. 0 i 0 


Los resultados de la estructura de bandas obtenidas por diagonalizar el hamiltoniano 
(4.36) se presentan en las figuras 28] y [29] para diferentes anchos del nanoribbons. La 
figura [28] corresponde a nanoribbons N-agnr, donde el número entero N que representa 


el número de dímeros que constituyen la cadena unidad del nanoribbons verifica que 


N+1 





M = es un número entero. (4.37) 


Los resultados muestran que los nanoribbons N-agnr donde N verifica la condición 
(4.37) tienen comportamiento semimetálico con gap de energía cero, debido a que en 
las proximidades de E = O tenemos al menos un estado permitido, sin importar el 


ancho del nanoribbons. 


Para los nanoribbons arnchair N-agnr, en el cual el número de dímeros N no cumple 
con la condición (4.37) la estructura de bandas se muestra en la figura [29] resultados 
en los cuales se observa que existe un gap de energía AE %X 0 sin importar el ancho de 


los nanoribbons, por tanto dichos nanoribbons tienen comportamiento semiconductor. 
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Figura 28: Espectro de energía de nanoribbons armchair cuyo ancho es formado por 


17-agnr, (e) 20-agnr y 


(d) 


agnr, (b) 11-agnr, (c) 14-agnr, 
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Figura 29: Espectro de energía de nanoribbons armchair cuyo ancho es formado por 
diferente cantidad de átomos (a) 9-agnr, (b) 10-agnr, (c) 12-agnr, (d) 13-agnr, (e) 18-agnr y 
(£) 100-agnr. 


4.3. Nanoribbons barbados (bearded) 


Para los nanoribbons de borde bearded mostrado en la figura eny=0y 
y = W, donde W = (3N — 1)a/2, donde N es el número de dímeros que consituyen 
la cadena unidad, encontramos que el hamiltoniano tight-binding se expresa por una 


matriz 2N x 2N obtenida de manera similar a la ecuación (4.18) como: 
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0 Í 0 0 0 0 
t 0 2t cos(k.a) 0 0 0 
0 2tcos(k,a) 0 de 0 0 
H=|0 0 Ñ 0 2tcos(k,a) 0 (4.38) 
0 0 0 2t cos(k.,a) 0 t 
0 0 0 0 l 0 


y 


lio liliolilta 














Figura 30: Espectro de energía de nanoribbons de borde zigzag-bearded con diferente 
número de átomos por celda unidad (a) 8-bearded, (b) 9-bearded, (c) 10-bearded, (d) 
11-bearded, (e) 30-bearded y (f) 100-bearded 
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La estructura de bandas de nanoribbons de borde beardeb obtenidos al diagonalizar 
el hamiltoniano (4.38) se presenta en la figura para diferentes anchos del 
nanoribbons, los resultados muestran que el espectro de energía es equivalente a los 
nanoribbons de borde zizag, pero difieren su comportamiento para energías próximas 
a E=0. 


4.4. Nanoribbons zizag-bearded 


Calculamos el espectro de energía correspondiente a nanoribbons con borde zigzag 
en y = W, y borde barbado (bearded) en y = 0, representamos por N-zgbearded, de 
modo que el número de átomos que constituyen el ancho W = 3NVa/2 del nanoribbons 
es 2N +1. 


La figura representa el espectro de energía para nanoribbons N-zgbearded, 
de diferentes anchos, obtenidas al diagonalizar el hamiltoniano de dimensión 
QN +1) x QN + 1). Los resultados encontrandos muestran que estos nanoribbons 
presentan un gap de energía AE 4 0 al rededor de E = 0 para cualquier ancho, es decir, 


los nanoribbons de borde mixto zigzag-bearded tienen comportamiento semiconductor. 
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Figura 31: Espectro de energía de nanoribbons de borde zigzag-bearded con diferente 
número de átomos por celda unidad (a) 7-zgbearded, (b) 8-zgbearded, (c) 9-zgbearded, (d) 
10-zgbearded, (e) 30-zgbearded y (f) 100-zgbearded 


4.5. Conductancia eléctrica 


De acuerdo a los espectros de energía observados en las figuras 
y los nanoribbons de grafeno presentan cuantización de estados permitidos [64], 
esta cuantización trasladada a la conductancia eléctrica produce la cuantización de 
la conductancia eléctrica que depende directamente del ancho del nanoribbons. Para 
el cálculo de la conductancia eléctrica utilizamos la formula de Landauer y las 
funciones de Green fuera de equilibrio [65]. A continuación presentamos resultados de 


conductancia eléctrica que afirman la naturaleza conductora o semiconductora de los 
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nanoribbons predicha por el espectro de energía. 


En la figura [32] mostramos los resultados obtenidos para la conductancia eléctrica, 
la densidad de estados y el espectro de energía de nanoribbons con borde zigzag, siendo 
la figura [32H(a) correspondiente a un nanoribbons 8-zgnr de ancho W = 11a, donde 
para energías —0.48t < E < 0.48t tenemos T(E) = 1, lo cual implica que para este 
rango de energía existe un solo canal (p, de transmisión permitido, esto se corrobora 
con el hecho de que para este rango de energía en el espectro de energía existe un 
solo valor de k permitido, en cambio para el rango de energía 0.48t < E < 0.76t 
tenemos T(E) = 3, que indica que, en este rango de energía energía existen tres 
canales (p, permitidos correspondientes a tres valores de k permitidos. El valor máximo 
del coeficiente de transmisión para el nanoribbons 8-zgnr es T(E) = 8, que ocurre 
en los rangos de energía 1.0t < E < 1.1t y —1.1t < E < 1.0t, asímismo en las 
energías E = +0.48t, +0.76t, +0.93t,...,+2.62t, +2.83t, +2.95t se observa una alta 


concentración de estados permitidos, que generan discontinuidades en la densidad de 








estados conocidos como discontinuidades de Van-Hove, y permiten el salto de 7 (E) en 
números enteros entre 0 y 8, ocurriendo el valor O para |E| > 2.95t. Así mismo en la 
figura [321(b) para el nanoribbons 12-zgnr de ancho W = 17a en el rango de energía 
—0.34t < E < 0.34t tenemos T(E) = 1, es decir un canal de transmisión, los valores 
permitidos para 7(E) en este caso van de O a 12, ocurriendo el valor O para |E| > 2.98t 
y el valor máximo 12 en los rangos de energía —1.05t < E < 1.00t y 1.00t < E < 1.05t. 
La naturaleza conductora de los nanoribbons con borde zigzag se debe al hecho de que 
para valores de energía próximas a E = O independiente del ancho del nanoribbons 
tenemos 7 (E) = 1, es decir T(E) es diferente de cero. 
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Figura 32: Conductancia, DOS y bandas de energía para (a) 8-zgnr y (b) 12-zgnr. 


Los coeficientes de transmisión obtenidas usando la fómula de Landauer, la densidad 
de estados y el espectro de energía de un nanoribbons de borde armchair 11-agnr 
de ancho W = 10y3a se presentan en la figura donde observamos que para 
energías —0.41t < E < 0.41t existe un solo canal p, permitido, ello debido a que 
T(E) = 1 y que en el espectro de energía k tiene un solo valor permitido. Para 
el rango de energía 0.41t < E < 0.48t existen dos canales (p, permitidos, debido 
a que T(E) = 2 y k tiene en este rango de energías dos valores perimitidos, el 
valor máximo del coeficiente de transmisión es T(E) = 5 que ocurre en el rango 
de energías —1.52t < E < —0.93t y 0.93t < E < 1.52t. Las discontinuidades de 
van-Hove que permiten el salto de T(E) en números enteros se presentan en las energías 
E = +0.41t; +0.48£; +0.72t; +0.93; +1.00€;...;+2.93t. En la figura [331b mostramos 


el coeficiente de transmisión, DOS y el espectro de energía correspondiente a un 








nanoribbons 14-agnr, que presenta propiedades similares al nanoribbons 11-agnr, con 
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la diferencia que los valores permitidos de T (E) varían de 0 a 7. 
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Figura 33: Conductancia, DOS y bandas de energía para (a) 11-agnr y (b) 14-agnr. 


Los resultados del coeficiente de transmisión y densidad de estados mostrados en 
la figura [33| verifican la naturaleza semimetálica de los nanoribbons armchair, cuyos 
anchos verifican la condición debido a que en las proximidades de E = 0 el 
coeficiente de transmisión es diferente de cero, es decir T(E) = 1 independientemente 


del ancho del nanoribbons. 


El coeficiente de transmisión, DOS y el espectro de energía de un nanoribbons de 
borde armchair 12-agnr de ancho W = 1143a/2 se presentan el la figura [3ala, donde 
observamos que en el rango de energías —0.13t < E < 0.13t tenemos T(E) = 0, en el 
rango de energías 0.13£ < E < 0.30t tenemos 7 (E) = 1, asímismo los valores de 7'(£) 
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son enteros que varían de O a 6. La figura [34+b muestra el coeficiente de transmisión 
) 


DOS y el espectro de energía de un nanoribbons 15-agnr, donde T7(E) varía de 0 a 7. 


(a) (b) 
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Figura 34: Conductancia, DOS y bandas de energía para (a) 12-agnr y (b) 15-agnr. 


El ancho de los nanoribbons armchair utilizados en la figura no cumplen con 
la condición estos nanoribbons tienen comportamiento semiconductor, debido a 
que T(E) = 0 en las proximidades de E = 0, es decir existe un pequeño gap de 
energía AE 4H 0. entre la banda de conducción y valencia. El gap de energía AF tiene 
una dependencia inversa con el ancho del nanoribbons, es decir a mayor ancho del 


nanoribbons el gap de energía disminuye. 


Para los nanoribbons barbados 8-bearded de ancho W = Sa los resultados del 


coeficiente de transmisión, DOS y el espectro de energía se muestran en la figura [35ta, 
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observandose que en el rango de energías —0.54t < E < 0.54t tenemos T(E) = 1, en el 
rango de energías 0.54t < E < 0.82t tenemos T(E) = 3, así mismo los valores de T'(£) 


varían de O a 7, ocurriendo el valor máximo en el rango de energías 0.98t < E < 1.18f£. 


(a) (b) 
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Figura 35: Conductancia, DOS y bandas de energía para (a) 8-bearded y (b) 12-bearded. 


Los resultados del coeficiente de transmisión, la DOS y el espectro de energía 
calculados para el nanoribbon de borde mixto zigzag-bearded 8-zgbearded de ancho 
W = 12a se presentan en la figura S6ha, donde se observa que en el rago de energía 
—0.38t < E < 0.38t tenemos T(E) = 0, existiendo por tanto un gap de energía 
AE = 0.76t entre la banda de condución y valencia. En la figura [36Fb para el 
nanoribbons 12-zgbearded encotramos que en el rango de energía —0.28£ < E < 0—28t 
tenemos 7 (E) = 0, así mismo se encuentra que 7'(E) varía entre O y 10, siendo el gap de 
energía entre banda de conducción y valencia AE = 0.56. Estos resultados corroboran 
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el hecho de que estos nanoribbons tienen comportamiento semiconductor, donde el gap 


de energía tiene una relación inversa con el ancho del nanoribbons. 
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Figura 36: Conductancia, DOS y bandas de energía para (a) 8-zgbearded y (b) 
12-zgbearded. 


4.6. Nanoribbons en campo magnético 


En esta sección describimos el espectro de energía de nanoribons de grafeno de 
diferentes bordes bajo el efecto de un campo magnético uniforme B perpendiccular al 
plano del nanoribbons, cuyo efecto es introducido en la integral de Hopping t;j por 


medio de la sustitución de Pierls 
do Bi (4.39) 
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donde a; es el flujo magnético Y por unidad de flujo cuántico Py = h/e dado por 


e f? 


Para un nanoribbons de borde zigzag N-zgnr la sustitución de Pierls (4.39) produce 
en los elementos de matriz 2t cos(k.,,a) de la ecuación (4.18) un cambio de fase de 27moa 


com m números enteros que van de 1 a NV, es decir 


2t cos(k.a) > 2tcos(k,a + 217ma), (4.41) 
donde 
== $ BoShex 
IT 


siendo Shez = 3/3 q? es el área de la celda unitaria de cada sub red del grafeno y By la 


intensidad del campo magnético aplicado. 


El espectro de energía de un nanoribbon 20-zgnr bajo la acción del campo magnético 
uniforme se muestra en la figura|[37|para diferentes valores de a = 9/0, (a) a = 0.0, (b) 
a; = 1/500, (e) a: = 1/250, (d) a = 1/100, (e) a = 1/50, (f) a. = 1/20, (8) a = 1/10, 
(h) a = 1/4 y (1) a = 1/2, encontrándose que para pequeños campos magnéticos 
(a = 1/500) el espectro de energía sufre ligeras variaciones respectro al espectro de 
energía para a, = 0, debido a que el movimiento ciclotrónico del electrón es interrumpido 
por el borde del nanoribbons, a medida que se incrementa la intensidad del campo 
magnético el movimiento ciclotrónico se encuentra mas definida, para a = 1/50 el 
radio ciclotrónico es comprarable con el ancho del nanoribbons. En cambio si el campo 
magnético es extremadamente alto a = 1/4 la estructura de bandas muestra que los 


niveles de Landau constituyen sub bandas de Landau. 
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Figura 37: Estructura de bandas de energía de un nanoribbons 20-zgnr (a) a: = 0.0, (b) 
a. =1/500, (c) a =1/250, (d) a = 1/100, (e) a = 1/50, (£) a = 1/20, (8) a = 1/10, (h) 
o=1/4y 0) a=1/2 


Para los nanoribbons de borde armchair 20-agnr y 22-agnr, los espectros de 
energía para diferentes valores de campo magnético se presentan en las figuras 
y respectivamente, donde similarmente a los nanoribbons 20-zgnr se observa que 
para pegueños campos magnéticos a = 1/500 los espectros de energía son casi 
invariantes respecto al espectro sin campo magnético a; = 0, resultado por tanto que 
los movimientos ciclotrónicos de las cargas eléctricas son interrumpidas por los bordes 
del nanoribbons, por tanto los niveles de Landau no se encuentran bien definidas, 


para campos magnéticos muy altos encontramos tambien que los niveles de Landau se 
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convierten en subandas de Landau. 


E/t 
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Figura 38: Estructura de bandas de energía de un nanoribbons 20-agnr (a) a = 0.0, (b) 
a = 1/500, (c) a = 1/250, (d) a: = 1/100, (e) a = 1/50, (f) a = 1/20, (8) a = 1/10, (h) 
u=147+(100a=1/2 
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Figura 39: Estructura de bandas de energía de un nanoribbons 22-agnr (a) a. = 0.0, (b) 
a: =1/500, (c) a = 1/250, (d) a = 1/100, (e) a = 1/50, (£) a = 1/20, (g) a = 1/10, (h) 
a=1Ay (M)a=1/2 


Cuando el campo magnético uniforme actúa sobre nanoribbons de borde bearded y 
bordes mixtos zigzag-bearded los espectros de energía para diferentes intensidades del 
campo magnético se presentan en las figuras [40] y [41] respectivamente. Los resultados 
muestran nuevamente que para pequeños campos magnéticos los movimientos 
ciclotrónicos de las cargas eléctricas se encuentran interrumpidas por los bordes de 
los nanoribbons, y para campos magnéticos extremadamente altos + = 1/4 los Niveles 


de Landau constituyen subbandas de Landau. 
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Figura 40: Estructura de bandas de energía de un nanoribbons 20-bearded 
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Figura 41: Estructura de bandas de energía de un nanoribbons 20-zgbearded (a) a. = 0.0, 
(b) a. = 1/500, (c) a = 1/250, (d) a = 1/100, (e) a = 1/50, (f) a =1/20, (g) a« = 1/10, (h) 
o=14y (0) 0 =1/2 
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Capítulo 5 


Nanoporos 


Los nanoporos son agujeros de tamaño nanométrico en materiales biológicos o 
de estado sólido que permiten transportar a través de él iones y moléculas, son 
dispositivos que han generado recientemente intensa actividad de investigación teórica 
y experimental en la comunidad científica. Los nanoporos se presentan en la actualidad 
como dispositivos muy prometedores para la exploración del genoma, en particular en 
la secuenciación del ADN [661167]. Se encuentra también que el principal inconveniente 
para utilizar nanoporos tanto en estado sólido o biológicos es la longitud (largo) es 
10 a 100 veces mayor que la distancia entre dos bases de una molécula de ADN (ver 
figura[42)a,b,c), esto hace que la lectura de los cambios de corriente eléctrica detectados 
por el nanoporo correspondan a un conjunto de moléculas de ADN y no así a una sola 
molécula, esta dificultad se salva si se utilizan nanoporos bidimensionales ultra delgados 
como es el caso de nanoporos de grafeno (ver figura [(42+d), que debido a su espesor 
de ún átomo permiten observar cambios característicos en la conductancia eléctrica 
producidos por el paso de una sola molécula a través del nanoporo. Para obtener los 
nanoporos de grafeno se colocan hojas de grafeno en agujeros de nitruro de silicio y se 
perfora el agujero de tamaño nanométrico en el grafeno utilizando un haz de electrones 
de un microscopio eléctrónico de transmisión [68]. Otra de las características que hacen 


interesante a los nanoporos de grafeno es que son materiales altamente conductores 


En los últimos años los nanoporos se han presentado como herramientas poderosas 
para el análisis de moléculas individuales, ellos han sido utilizados para caracterizar 
rápida y eficientemente biopolímeros de ADN usando nanoporos de estado sólido 
[69)[70], se examinaron también el efecto local de proteínas a lo largo de la longitud 
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del ADN, usando el hecho de que el transporte de corriente a través el ADN, así 
como del ADN con proteínas a lo largo de su longitud son distintas [71]. Estos y 
otros resultados experimentales son un paso crucial para la exploración del genoma, y 
permiten visualizar la viabilidad de la lectura de la información a lo largo del ADN a 


alta resolución con un nanoporo de estado sólido es posible. 
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Figura 42: Nanoporos para la secuenciación del ADN. (a) el ADN se inserta en un 
nanoporo, con control de velocidad proporcionado por una ADN polimerasa 929 (marrón). 
Un nanoporo a-hemolisina (gris) está incrustado en una bicapa de lípidos. El sustrato de 
ADN (columna vertebral rojo) se inserta en el poro por un campo eléctrico aplicado. Su 
movimiento dentro o fuera del poro (flechas) puede ser controlado por el campo eléctrico 
aplicado y la actividad de la polimerasa. (b) Estructura de un nanoporo a-hemolisina 
insertado en una bicapa lipídica. (c) Estructura de un nanoporo MspA insertado en una 
bicapa lipídica. (d) Un nanoporo grafeno dibujado en la misma escala que la a-hemolisina y 
poros proteicos MspA [72]. 


En la actualidad encontramos también intensa actividad científica en el trasporte 
electrónico de sistemas moleculares principalmente en el caso del ADN, por su 
importancia médica y la posibilidad de construir nuevos dispositivos que utilicen 
sistemas híbridos de ADN con materiales de estado sólido, debido a que los 


resultados encontrados muestran que las moléculas de ADN tienen características 
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semiconductoras, que depende del segmento de ADN y en particular de la secuencia 
de las moléculas de ADN [73476]. Estos resultados permiten también entender que la 
utilización de nanoporos de grafeno es perfectamente posible para la lectura directa de 
diversos tipos de moléculas por ser el grafeno un material conductor por excelencia [77]. 
Los cambios que generan cada tipo de las secuencia nitrogenada de ADN, conocidos 
también como nucleótidos, que son la adenina, guanina, citosina y timina producen 
al pasar por los nanoporos un cambio característico en la conductancia eléctrica del 
nanoporo, cualidad que permitirá en el futuro realizar una lectura rápida y eficaz de 


secuencias completas del ADN. 


Nuestra idea en el presente capítulo es analizar numéricamente dentro del modelo 
del tight-binding y el régimen de transporte balístico los cambios que producen diversos 
tipos de moléculas en los perfiles de conductancia eléctrica de nanoporos construidos 
en redes cuadradas y el grafeno, la descripción del transporte electrónico se realiza 
usando la Formula de Landauer y las funciones de Green descritas en el capítulo 3. 
Los resultados se muestran en términos de la variación del coeficiente de transmisión, 
debido a que los cambios producidos en los perfiles del coeficiente de transmisión por 


la presencia de las moléculas son muy pequeñas. 


5.1. Nanoporos en red cuadrada 


El transporte de carga a través de nanodispositivos adheridos a dos electrodos 
metálicos semi-infinitos se describe en la sección donde la Fórmula de Landauer 
establece que la conductancia eléctrica se encuentra cuantizada en múltiplos de la 


unidad cuántica de la conductancia Gy = e*/h. 


En los últimos años, el transporte de carga a través de sistemas moleculares 
ha despertado el interés de diferentes grupos de investigación. Y obviamente las 
moléculas que despiertan mayor interés son los cuatro monómeros o nucleótidos del 
ácido desoxirribonucleico, ADN, que consisten de dos partes, la primera es un azúcar, 
desoxirribosa, con un grupo fosfato adherido, mientras que la segunda es una base que 
puede ser la adenina (A), guanina (G), citosina (C) o timina (T) [78]. El esqueleto 
del ADN está formado por dos cadenas que forman los azúcares unidos a través de los 
grupos fosfatos y las bases se colocan perpendiculares a las cadenas de modo que unen 
a las mismas pero sin repetir la base en el mismo lugar de las dos cadenas, las bases se 
unen como CG y AT, o GC y TA a través de dos pares ligados de hidrógeno [78]. 


TI 


Se han reportado estudios teóricos sobre el transporte de carga a través de algunas 
secuencias de ADN que pueden ser considerados desde el punto de vista de la 
física como nanohilos conectados a electrodos que en algunos estudios son nanoribbons 
de grafeno [75]. Los resultados sustentan que esas secuencias de ADN tienen un 
comportamiento semiconductor y muestran que la relación corriente-voltage de los 


segmentos depende de la longitud del nanohilo y en particular de la secuencia de los 


nucleótidos , 


Por otro lado, se ha sugerido que cuando un sistema conductor no son moléculas 
sino una red cristalina bidimensional como el grafeno con un nanoporo conectado a dos 
contactos métalicos semi-infinitos debe presentar una conductancia que podría sensar 
moléculas que atraviesen el nanoporo de modo que la conductancia dependiese del tipo 
de molécula y la orientación de la misma al atravezar el nanoporo. De esa forma se 
podría identificar los nucleótidos del ADN sensando las bases que presentarían una 


seña particular y característica en la conductancia para cada base [77]. 


En la presente sección, hacemos un estudio bastante básico y fundamental de un 
nanoporo en un conductor rectangular adherido a dos contactos metálicos semi-infinitos 
donde el conductor y los contactos son emulados por una red cuadrada en la 
aproximación tight-binding a primeros vecinos. Sobre el nanoporo colocamos moléculas 
con orientaciones en el eje x e y y consideramos que las moléculas están débilmente 


acoplados al conductor. 


5.1.1. Modelo 


Un conductor adherido a dos contactos metálicos semi-infinitos con potenciales 
químicos determinados es un sistema abierto. Esto significa que los portadores de 
carga que proceden de izquierda a derecha con una energía igual al potencial químico 
del contacto izquierdo son esparcidos por las interfaces contacto izquierdo-conductor y 
conductor-contacto derecho, de modo que al alcanzar el contacto derecho, los portadores 
de carga no regresan hacia el contacto izquierdo. De la misma forma, los portadores de 
carga que proceden de derecha a izquierda -con una energía igual al potencial químico 
del contacto derecho-, sufren esparcimiento en las interfaces mencionadas y una vez 


alcanzado el contacto izquierdo no retornan más hacia el contacto derecho. 


Considerando el punto de vista reduccionista, los contactos y el conductor pueden 
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ser emulados por redes rectangulares semi-infinitas de sítios atómicos en el primer 
caso y por una red rectangular finita de sítios atómicos en el segundo. Esto es posible 
porque, cuando un portador de carga se mueve en una dimensión el hamiltoniano 
está conformado por la energía cinética y la energía potencial que sufre el portador. 
Cuando se discretiza la energía cinética se genera una matrix tridiagonal y cuando se 
hace lo mismo con la energía potencial se obtiene una matriz diagonal. Los elementos 
no diagonales de la matriz de energía cinética está dado por 
h? 

e (5.1) 
donde m* es la masa efectiva de los portadores de carga, a es el parámetro de red 
de la cadena lineal y Ah la constante de Planck dividido entre 27. Y los elementos 
diagonales de la matriz están dados por U; + 2t, donde U; es la energía potencial del 
punto de discretización j y a t se le denomina el término de acoplamiento entre los 
puntos que resultan de la discretización de un sistema unidimensional (1D) [46]. Esta 
matriz producida por la discretización por diferencias finitas del hamiltoniano de un 
flujo de portadores de carga a través de un conductor unidimensional es equivalente 
al hamiltoniano de una cadena unidimensional de sitios atómicos que se acoplan 
considerando orbitales atómicos tipo s a primeros vecinos en la aproximación de 


ligaciones fuertes (tight-binding) con un parámetro de hopping t [46]. 
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Figura 43: Representación esquemática de un nanoporo rectangular construído en una red 


rectangular finita conectado a dos electrodos semi-infinitos. 


Por otro lado, para una red rectangular finita de sítios atómicos con orbitales tipo 
s y acoplamiento a primeros vecinos, el hamiltoniano en primera cuantización está 


definido como 
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1,3 UN) 


donde t.. y t,y son los parámetros de hopping a lo largo del eje x e y, respectivamente, 


y los índices 1,7 indican las posiciones discretas del los sítios de la red dadas por 


(e y;) — (10. Oy): 


Un nanoporo es generado en la red rectangular finita considerando la energía de los 
sítios en esa región un valor totalmente diferente y mayor que la energía de los sítios 
de los contactos izquierdo y derecho, pz y pp, respectivamente y de la red rectangular 
HL — HR, como muestra el diagrama esquemático de la figura donde el nanoporo 
de longitud D y ancho W es construido dentro de la red finita de longitud X y ancho 
Y, que está conectada a los electrodos izquierdo y derecho que se identifican con los 
potenciales químicos 4, y pr, respectivamente, el electrodo izquierdo constituye la 
fuente de cargas eléctricas y el electrodo derecho el sumidero. El análisis del transporte 
de portadores se estudia dentro de la formulación de Landauer descrito en la 
sección donde la conductancia eléctrica G se expresa en función del coeficiente de 
transmisión 7' (E) como 


au)=*ND) (5.3) 


En este contexto, tenemos el primer proceso en que definimos la conductancia del 
sistema contacto izquierdo-conductor-contacto derecho; luego el segundo proceso en la 
que obtenemos la conductancia del sistema con un nanoporo incluído en el conductor 
y el tercer proceso, que involucra la presencia de moléculas en el centro del nanoporo 


y calculamos la conductancia del sistema. 


Para simular el efecto de las moléculas -que se encuentran en el centro del nanoporo-, 
sobre la conductancia del sistema construímos un algoritmo que incluye el hamiltoniano 
de la molécula en el nanoporo considerando la interacción entre átomos de la propia 
molécula, así como la interacción entre la molécula con los sítios del borde que definen 


el nanoporo dentro de la aproximación tight-binding. 


80 


Funciones de Green 


Para determinar los canales de transmisión en el dispositivo resolvemos la ecuación 
de autovalor del Hamiltoniano del sistema usando el método de las funciones de 
Green [161/65). Finalmente, la función de Green Gp del conductor con el nanoporo 
y las moléculas en el centro del nanoporo tal como se muestra en la figura [43] se obtiene 
como 


i 
[E bs re (5.4) 


donde / es la matriz identidad, Hp el hamiltoniano del conductor que incluye al 
nanoporo y a las moléculas que se incluyen en el centro del mismo, +in > 0 es un 
número imaginario pequeño que permite obtener las funciones de Green retrasadas 
(-) y avanzadas (+); asimismo 2, son las auto-energías asociadas a los electrodos 
izquierdo y derecho, respectivamente, y que contienen toda la información acerca de 


los contactos semi-infinitos. 


Para determinar las auto-energías de los contactos semi-infinitos usamos las ideas 
propuestas por Lee y Joanmnopoulos quienes consideran que los conductores 
semi-infinitos están formados por cadenas verticales finitas de sítios, hy, que se acoplan 
a las cadenas verticales laterales y vecinas a izquierda y a derecha mediante una 
matriz diagonal, V, ,, cuyos elementos incorporan las interacciones a primeros vecinos 
entre cada sítio de las cadenas verticales. Por tanto, las funciones de Green de los 
contactos izquierdo y derecho están dados por q, = Vi(1/A)V* y q, = V¿(1/h,)V7, 


respectivamente. 


Luego las auto-energías de los contactos se determinan mediante 2%, = 2, siendo 
a = Hny94H ¡ siendo H py la matriz de acoplamiento entre el conductor y el contacto 
q y Ha = Hop la matriz de acoplamiento del contacto q con el dispositivo. Para el 


cálculo de gr utilizamos el método recursivo de convergencia rápida sugerida por 


Sancho et al [50151]. 


Finalmente, el coeficiente de transmisión es 7(E) = Tr P¿GpPrGÍ| donde Pz r 
son las funciones de acoplamiento de los contactos con el dispositivo y se determinan 


mediante Pp = (Lor — Bla) 


A partir de la función de Green, Gp, obtenemos también, la densidad de estados 
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(DOS) que proporciona la cantidad de estados disponibles en un intervalo de energía, 


p(E) = -m [In [Gp(2)||. (5.5) 


5.1.2. Resultados y discusión 


En el presente trabajo simulamos el conductor como una red rectangular finita de 
sítios que poseen una periodicidad en las dos dimensiones de modo que los parámetros 
de red son 4, = y considerando una interacción a primeros vecinos con una energía 
de acoplamiento t, = t,. Los contactos semi-infinitos conectados a ambos lados del 
conductor poseen la misma periodicidad espacial y la misma interacción energética a 
primeros vecinos. Debido a las reducidas dimensiones del sistema, los coeficientes de 
transmisión a través del conductor se cuantizan de forma tal como se observan en la 


figura 


La figura[44[a) muestra la evolución de la cuantización del coeficiente de transmisión 
con respecto a la energía y se puede apreciar que los valores del coeficiente de tranmisión 
evolucionan desde O hasta 20. El valor O corresponde para los coeficientes de transmisión 
cuando la energía, |E| > At, y- Para estos valores de energía no existe transporte porque 
la energía de los portadores de carga son superiores o inferiores al ancho de la banda del 
sistema. Esto se corrobora observando las densidad de estados, DOS, cuyos resultados 
se muestran en la figura ), observándose que para |E| > 4£,, y no existe presencia 


de estados permitidos que permitan el transporte de portadores de carga. 


El coeficiente de transmisión máximo de 20 ocurre cuando la energía de los 
portadores alcanza el nivel del potencial químico E = Ep = 0. Asimismo, observamos 
que los saltos en la cuantización del coeficiente de transmisión ocurren en valores 
de energía donde la DOS presenta las discontinuidades de Van-Hove. Los resultados 
encontrados para el coeficiente de transmisión y la DOS de la red finita son 
equivalentes con los coeficientes de transmisión y la DOS de una nanocinta formada 
por una red cuadrada de ancho y = 20a,,, debido esencialmente a las condiciones de 
transporte balístico asumidas en el formalismo de Landauer, en la que los procesos 
de esparcimiento (scattering) de los electrones en la red rectangular finita ocurren 
únicamente en las interfaces entre la red finita y los electrodos, dependiendo la 


cuantización del coeficiente de transmisión solamente del ancho y. 


Los efectos de la presencia de los nanoporos en el conductor rectangular en el 
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transporte eléctronico se muestran en las figuras [44[c), (e) y (g), donde se observan 
que el coeficiente de transmisión disminuye con respecto al coeficiente de transmisión sin 
poros Fig. [44(a), dicha disminución se acentúa a medida que el ancho W del nanoporo 
aumenta. También se observa que la cuantización en números enteros de e*/h del 
coeficiente de transmisión se va perdiendo, además se observa que el coeficiente de 
transmisión es independientes de la longitud D del nanoporo, corroborando el aspecto 


balístico del transporte. 





Energía / t Energía / t 


Figura 44: Se muestran los coeficientes de transmisión (izquierda) y DOS (derecha) para una 
red rectangular sin poros (a) y (b); para redes rectangulares con poros de D=9a,, y W=3ay , 
(c) y (d); D=94. y W=5ay, (e) y (£); y D=907 y W=7ay, (8) y (h), respectivamente. 
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Los resultados para las DOS que muestran las figuras Lala), (£) y (h) indican una 
disminución progresiva respecto a la DOS sin poros, figura a), principalmente en 
aquellos valores de energía donde la DOS presenta las discontinuidades de Van-Hove, 
es decir, que el número de estados permitidos en el sistema disminuye por la presencia 


del nanoporo y el aumento del tamaño del nanoporo en la dirección y. 


A continuación escogemos un nanoporo cuadrado de dimensiones D=W=9a,,, en el 
centro del mismo simulamos la presencia de una molécula unidimensional orientadas 
en la dirección x y luego otra molécula en la dirección y tal como se observa 
esquemáticamente en la Figura Consideramos que las moléculas presentan una 
interacción entre sus átomos en el esquema tight-binding con parámetro de hopping t.,, 
de modo que el hamiltoniano molecular está dado por H,, = > tm[10)(6+1|+]/|¿+1) (1), 
siendo |¿) los autoestados de energía asociados a los orbitales tipo s de los sitios de los 
átomos. Asimismo la interacción entre el nanoporo y la molécula se asume, también, 
dentro del modelo tight-binding, de modo que un átomo de la molécula interactúa con 


el primer sitio más próximo del nanoporo por medio del parámetro de hopping tmp 








dado por 
1? h? te 
tmp =- =- = (5.6) 
2mb? 2M(Tzyle y)? Toy 


siendo hb = 1, y; y la distancia en la dirección paralela x y perpendicular y a la corriente 


entre la molécula y el nanoporo. 
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Figura 45: Representación esquemática de un nanoporo cuadrado con D=W y una molécula 
de cuatro átomos colineales orientados en dirección paralela (a) y en la dirección perpendicular 


(b) al flujo de la corriente. 


Los resultados para el coeficiente de transmisión del sistema con un nanoporo de 
dimensiones D=W=09a,, en la cual se encuentra una molécula lineal de cuatro átomos 
con acoplamiento t,, = t se muestran en la figura donde las curvas en negrita 


representan los coeficiente de transmisión del sistema con nanoporo y la presencia de la 
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molécula, mientras que las curvas azules corresponden a los coeficientes de transmisión 
del sistema con nanoporo sin la presencia de moléculas. La figura 46[a) corresponde al 
coeficiente de transmisión asociado al nanoporo con una molécula orientada paralelo 
al flujo de la corriente, Figura [45(a), donde la presencia de la molécula en el nanoporo 
produce resonancias (ver en el ¿nset de la Figl46[a)) en el coeficiente de transmisión 
en cuatro valores de energía E/t = +0.62; +1.62, que corresponden al acoplamiento 
entre los átomos que componen la molécula. En la figura b) observamos que el 
coeficiente de transmisión para el sistema con nanoporo y con molécula orientada 
perpendicularmente al flujo de corriente, Figura ), los coeficientes de transmisión 
próximos al potencial químico (T=0) o energía de Fermi presentan antiresonancias (ver 
inset de la Fig[46[b)) mientras que para energías no muy próximas al nivel de Fermi, 


las resonancias son muy bajas. 
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Figura 46: Se muestran los coefientes de transmisión en un conductor con un nanoporo 
rectangular de dimensiones D=W=09az, en la presencia de una molécula lineal con cuatro 
átomos orientada paralela (a) y perpendicularmente (b) al flujo de corriente eléctrica. En el 


primer caso se presentan resonancias y en el segundo anti-resonancias (ver insets). 


Debido a que las resonancias y antiresonancias que se observan en los coeficientes 
de transmisión son pequeñas, en las figuras siguientes representamos las variaciones del 
coeficiente de transmisión AT' que mide la diferencia del coeficiente de transmisión del 
sistema con nanoporo con la presencia de la molécula y del coeficiente de transmisión del 
sistema con nanoporo sin la molécula. Los resultados que se muestran en la Figura [47] 
corresponden a AT(E) para el sistema con un nanoporo con moléculas orientadas 
paralelas a la corriente de portadores de carga, de modo que las moléculas varían 
en el número de sus componentes atómicos. Los resultados para una molécula con 


tres átomos se muestran en la Fig. a) presentando tres picos en E/t = 0.0; +1.43; 
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para una molécula compuesta con cuatro átomos, la Fig. [47/b), muestra los picos 
correspondientes a cuatro energías E /t = +0.62; 1.62; para una molécula con 5 átomos, 
la Fig. c), presenta los picos en E/t = 0.0;+1.0;+1.73; y finalmente, para una 
molécula compuesta por seis átomos, la Fig. (47(d), presenta los picos en los seis valores 
de energías dados por E /t = +0.46; +1.25; +1.80. Estos resultados indican exactamente 


la distribución espectral de la energía de acoplamiento de un número dado de átomos 





que conforman la molécula, asimismo observamos que la altura de los picos son mayores 
a medida que el número de átomos de la molécula aumenta. 
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Figura 47: Variación del coeficiente de transmisión, AT, para un conductor con un nanoporo 
de dimensiones D=W=9a y con la presencia de moléculas lineales con varios átomos idénticos 


(a) 3 átomos, (b) 4 átomos, (c) 5 átomos y (d) 6 átomos, paralelos a la corriente de portadores. 


A diferencia de la figura [47] donde las moléculas lineales son paralelas al flujo de la 
corriente, en la figura [48| mostramos resultados correspondientes a moléculas lineales 
orientados en dirección perpendicular a la corriente de portadores, observandose que las 
resonancias y antiresonancias ocurren en cada caso para los mismos valores de energía, 


predominando en todos los casos las antiresonancias respecto a las resonancias. 
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Figura 48: Variación del coeficiente de transmisión, AT, para un conductor con un nanoporo 
de dimensiones D=W=9a y con la presencia de moléculas lineales con varios átomos idénticos 
(a) 3 átomos, (b) 4 átomos, (c) 5 átomos y (d) 6 átomos, perpendiculares a la corriente de 


portadores. 


Con el fin de emular el efecto en el nanoporo de moléculas diferentes a las moléculas 
lineales consideramos moléculas mostradas en la figura que se encuentran en el 
centro de un nanoporo rectangular de dimensiones D = W = 9a,, y teniendo en cuenta 
que cada sitio de la molécula interactúa con su vecino próximo con el parámetro de 
Hopping t,, = t, los espectros correspondientes para AT' se muestran en la figura [50] en 
las cuales los picos positivos indican valores de energía en las cuales la molécula produce 
un incremento en el transporte electrónico a través del nanoporo y los picos negativos, 
indican que en dichos valores de energía produce un aumento en la resistividad eléctrica, 


que disminuye el transporte electronico. 


(a) (b) (c) (d) 


Figura 49: Diferentes tipos de moléculas que se consideran en el centro del nanoporo 
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Figura 50: Variación del coeficiente de transmisión AT de nanoporo con D =W = Ya 


bajo la presencia de moléculas representadas en la figura [49] respectivamente. 


Conclusiones 


En un nanoporo rectangular construido en una red rectangular que emula al GaAs 
el transporte de carga depende del ancho del nanoporo, afirmación que se sustenta en 
el hecho de que el coeficiente de transmisión disminuye a medida que incrementa el 


ancho del nanoporo. 


Las moléculas presentes en el centro del nanoporo rectangular producen picos en 
el espectro de AT' que son característicos para cada tipo de molécula, dichos picos 
dependen del tipo de átomo, y número de átomos que constituyen la molécula, así 
mismo dependen de la forma en el cual se distribuyen los átomos en la molécula. 
Estas características proporcionan la posibilidad de utilizar los nanoporos de carbono 
como una técnica espectroscópica para el análisis de materiales, para tal fin se 
requiere construir una base de datos completa de cambios en la conductancia eléctrica 
(coeficiente de transmisión) producida por diferentes moléculas y desarrollar tecnicas 


experimentales que permitan medir los pequeños cambios en el transporte de carga. 
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5.2. Nanoporos de grafeno 


En esta sección analizamos el transporte electrónico en nanoporos rectangulares 
construidos en una red de grafeno considerando el efecto producido en el transporte 
electrónico debido a la presencia de moléculas con diferentes configuraciones en el centro 
del nanoporo, el comportamiento de los electrones es descrito utilizando el método los 
enlaces fuertes (tight-binding), el transporte electrónico se realiza en régimen balístico 


por medio de la formula de Landauer, y las funciones de Green. 


Como punto de partida para el estudio consideramos que el nanoporo de grafeno 
es un nanodispositivo conectado a dos electrodos metálicos semi-infinitos, que debido 
a sus pequeñas dimensiones tiene una conductancia eléctrica cuantizada en múltiplos 


de la unidad de conductancia cuántica Gp = e?/h. 


En los últimos años el transporte de carga a través de nanodispositivos moleculares 
ha despertado gran interés en diferentes grupos de investigación. Particularmente las 
moléculas que despiertan mayor interés son los cuatro monómeros o nucleótidos del 
ácido desoxirribonucleico, ADN, que consisten de dos partes, la primera es un azúcar, 
desoxirribosa, con un grupo fosfato adherido, mientras que la segunda es una base que 
puede ser la adenina (A), guanina (G), citosina (C) o timina (T) [78]. El esqueleto 
del ADN está formado por dos cadenas que forman los azúcares unidos a través de los 
grupos fosfatos y las bases se colocan perpendiculares a las cadenas de modo que unen 
a las mismas pero sin repetir la base en el mismo lugar de las dos cadenas, las bases se 
unen como CG y AT, o GC y TA a través de dos pares ligados de hidrógeno [78]. 
Como mencionamos al inició del capítulo los nanoporos de grafeno surgen surgen 
como una alternativa poderosa para el análisis de moléculas y el secuenciamiento del 
ADN, esto es posible gracias a que el espesor del grafeno es de un átomo, propiedad 
que permite leer pequeñas variaciones que sufre la corriente del nanoporo cuando 
una molécula pasa a través de nanoporo, estas variaciones son caracteristicos del 
tipo de molécula y sus átomos constituyentes. La secuenciación del ADN usando 
nanoporos ofrece grandes ventajas respectos a los métodos bioquímicos por ser un 
proceso de lectura de moléculas rápida y directa, y no requiere complejos procesos 
bioquímicos de crecimiento, clonado y reacciones químicas necesarios en los procesos 
bioquímicos. Varios grupos de investigadores an reportado muchas características del 
transporte de carga de moléculas de ADN, se conoce que cada nucleótido adenina 
(A), guanina (G), citosina (C) y timina (T') producen una relación corriente - voltage 
característico [681[72]177][81), los resultados muestran el trasporte de carga a lo largo 
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del esqueleto de segmentos de ADN y el transporte de carga en cada nucleótido en 
la dirección transversal al esqueleto del ADN. La finalidad de la secuenciación del 
ADN es la determinación del orden de los nucleótidos (A, C, G y T) que constituye 
la información genética heredable del núcleo celular en plantas y animales, por tanto 
es útil en el estudio de los procesos biológicos fundamentales que permiten determinar 
mutaciones, presencia de células cancerígenas, etc. Las primeras pruebas experimentales 


de translocación de moléculas de ADN fueron reportados en 1996 usando nanoporos 


biológicos [82]. 


La idea de utilizar nanoporos de grafeno en el proceso de análisis de moléculas 
individuales y la secuenciación del ADN es medir los perfiles de conductancia del 
nanoporo y observar las pequeñas variaciones que sufre el perfil del conductancia 
cuando una molécula o un nucleótido de ADN se encuentra en el centro del nanoporo, 
dichas variaciones del perfil de conductancia son característicos del tipo del molécula y 
permiten por tanto una lectura rápida y directa del tipo de molécula que se encuentra en 
el centro del nanoporo. En la actualidad se viene trabajando fuertemente en el desarrollo 
de técnicas de medición de corriente que sean sensibles a las pequeños cambios de 
corriente originadas por las moléculas y que permitan la producción comercial del los 


dispositivos de secuenciamiento de ADN, 


En el presente estudio la interacción de los electrones con los átomos del 
nanodispositivos se describe dentro del modelo de enlaces fuertes (método del 
tight-binding) a primeros vecinos y el transporte de carga se analiza en el régimen 
balístico, caso en el que la dispersión de las cargas eléctricas que ingresan al 
nanodispositivo desde los electrodos se realiza solamente en las superficies de contacto 


entre el nanodispositivo y los electrodos y en los bordes de los nanoporos. 


5.2.1. El modelo 


El diagrama esquemático de un nanoporo rectangular de grafeno conectado a 
electrodos ideales semi-infinitos con potenciales químicos 4, y pp se muestra en la 
figura donde consideramos que el electrodo izquierdo de potencial químico uy es 
la fuente de portadores de carga y el electrodo derecho con potencial químico ya es el 
sumidero. El transporte electrónico se considera balístico, es decir que los portadores 
de carga que proceden de los contactos hacia el nanoporo sufren dispersión solo en las 


interfaces entre el nanoporo y los contactos y en el borde del nanoporo, de modo que 
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toda carga eléctrica que procede del contacto izquierdo con energía igual a 1, y que 
alcanza el contacto derecho no regresa hacia el contacto izquierdo, de la misma forma 
toda carga eléctrica que ingresa por el contacto derecho con energía ug y alcanza al 
contacto izquierdo no volverá hacia el contacto derecho. 


contacto contacto 
izquierdo nanoporo derecho 





Hr «o X Pr 


Figura 51: Representación esquemática del modelo de nanoporo rectangular en grafeno 


conectado a dos electrodos semi-infinitos de borde zigzag 


De acuerdo a la figura el nanoporo rectangular de longitud D y ancho W es 
construido en una red finita de grafeno de longitud X y ancho Y, y se encuentra 
conectada a dos electrodos ideales a potenciales químicos uu, y Ur. El análisis del 
transporte de carga se efectúa dentro de la formulación de Landauer donde 
la conductancia eléctrica G se determina en función del coeficiente de transmisión 7 (4) 


como 


De? 
G= ÁT(Er), (5.7) 


Contexto en el cual nuestro primer paso es determinar la conductancia del sistema 
formado por la red finita de dimensiones X x Y de grafeno con dos contactos, el segundo 


paso es obtener la conductancia del nanoporo de dimensiones D x W construido en 
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la red finita de grafeno, y el tercer paso es calcular la conductancia considerando la 


presencia de moléculas de diferentes configuraciones en el centro del nanoporo. 


La interacción entre sitios atómicos en el nanoporo se consideran dentro de la 
interacción tight-binding, descrita en la sección 2.1 en términos de los estados cuánticos 


lo?) y |p5) asociados a las sub redes triangulares A y B del grafeno como 


H=-t) lp] +h.c., (5.8) 
(8,3) 
donde t = —2.8 eV es el parámetro de hopping y (i,j) indica que la suma se efectúa 


sobre los vecinos próximos. 


Para simular el efecto de las moléculas - que se encuentran en el centro del nanoporo 
- construimos un algoritmo que incluye en el hamiltoniano del nanoporo al hamiltoniano 
correspondiente a la molécula y la interacción de la molécula con el nanoporo. La 
interacción entre los átomos de la propia molécula, así como la interacción de la 
molécula con el borde que define al nanoporo se efectúa dentro de la aproximación 


tight-binding. 


El algoritmo que permite calcular la conductancia G del nanoporo utilizando la 
ecuación calcula primero las funciones de Green g, de los contactos izquierdo 
(q = L) y derecho (q = R), teniendo en cuenta que cada contacto se encuentra formada 
por cadenas verticales periódicas, donde cada cadena interactúa solo con las cadenas 
paralelas y vecinas a izquierda y derecha, esta propiedad permite obtener g, utilizando 
el método de convergencia rápida [50/51]. Segundo calcula las auto-energías Xy y las 
matrices de acoplamiento I',y = (2, — El) de los contactos. Tercero calcula la función 
de Green del sistema nanodispositivo y contactos usando las funciones de Green fuera 


de equilibrio, que de acuerdo a la sección [3.2.2] se expresa como: 


G=|[EI- Hp-Y,-—Yn]* (5.9) 


aquí / es la matriz identidad, Hp el hamiltoniano del nanodispositivo en estudio y 2, 
2 son las autoenergías asociadas a los contactos izquierdo y derecho respectivamente. 


Cuarto calcula el coeficiente de transmisión del sistema completo como 


T(E) =Tr[T¿GpTrG5), (5.10) 
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y la densidad de estados, DOS, que proporciona la cantidad de estados disponibles por 


cada energía permitida E como 
1 
p(E) = - TrlIm(G(E))]. (5.11) 


5.2.2. Resultados y discusión 


Simulamos primero una red finita de grafeno de longitud X = 20y3a y ancho 
Y = 28a con dos contactos, considerando interacción a primeros vecinos con parámetro 
de acoplamiento de Hopping t = —2.8 eV. Los contactos semi-infinitos se consideran 
cintas de grafeno de ancho Y con el mismo parámetro de Hopping que la red finita. 
Debido a las pequeñas dimensiones de la red finita los coeficientes de transmisión 
encontrados se encuentran cuantizados en valores enteros que van de O hasta 18 como 
se muestran en la figura |[52(a), la densidad de estados correspondiente a este sistema 
muestra que los saltos en el coeficiente de transmisión ocurren en las discontinuidades 
de Van-Hove como se observa en la figura B2[b). Los resultados del coeficiente de 
transmisión y la DOS son equivalentes al de un nanoribbons de borde zigzag 18-zgnr, 
debido a las condiciones de transporte balístico asumidas en el formalismo de 
Landauer, es decir el resultado de T(E) es independiente de la longitud X del sistema, 
dependiendo la cuantización solamente del ancho Y. El nanoribbons tiene naturaleza 
conductora debido a que en las proximidades de Ep = 0 tenemos que T(E) 4 1, es 


decir no existe presencia de un gap de energía entre la banda de conducción y valencia. 


Con la finalidad de observar los cambios que sufre el coeficiente de transmisión 
de la red finita de grafeno de dimensiones X x Y bajo la presencia de nanoporos, 
las dimensiones del nanoporo se muestran en la figura donde D es la longitud 
del nanoporo y W es el ancho del nanoporo. Los efectos producidos por la presencia 
de nanoporos rectangulares de diferentes dimensiones D x W sobre el coeficiente de 
transmisión 7 (E) se muestran en las figurasB52[c), (e) y (8), casos en las que se consideró 
D = 7vy3a y distintos anchos W = 4a, W = 7a y W = 10a respectivamente. El 
resultado muestra que el coeficiente de transmisión del sistema con presencia de poros 
es menor que el coeficiente de transmisión sin presencia de poros, dada en la figura 
B2la), esta disminución se acentúa a medida que el ancho W del poro aumenta, así 
mismo debido a las características de transporte balístico los resultados son invariables 


respecto a la longitud D del nanoporo, se observa tambien que la cuantización en 
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múltiplos enteros de la conductancia cuántica e*/h se va perdiendo. 
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Figura 52: Coeficientes de transmisión (izquierda) y Densidad de estados (derecha) para 


una red finita de grafeno con dos contactos y sin poros (a) y (b), para la red de grafeno con 


poros (c) y (d) D= 743a y W = 4a (e) y (f) D = 7V3a y W = Ta (8) 


W 


y (h) D=743a y 
= 10a 


Los resultados mostrados en las figuras [52[d), (£), y (h) muestran que la DOS del 


sistema con presencia de poros es menor que la DOS del sistema sin presencia de 
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poros, dada en la figura [52(a), la disminución se acentúa principalmente en aquellos 
valores de energía donde sin la presencia de poros se presentaban las discontinuidades 
de Van-Hove, esto se debe a medida que se incrementa el tamaño del poro el número de 
estados permitidos en el sistema disminuye. En los nanoporos rectangulares estudiados 
la disminución en la DOS del sistema con poros depende de la medida del ancho 
W del nanoporo y es independiente de la longitud D del nanporo, ello debido a las 


características de transporte balístico considerado. 





Figura 53: Representación esquemática de un nanoporo rectangular de grafeno con una 


molécula compuesta por cuatro átomos colineales 


A continuación dentro de nanoporo de longitud D = 743a y ancho W = 7a, cuyos 
coeficientes de transmisión y DOS son dados en la figura [44[c) y (d) respectivamente 
emulamos la presencia de una molécula lineal orientada en la dirección del eje z, 
esta molécula lineal es una red unidimensional finita como muestra la figura 
Consideramos que la interacción entre sitios atómicos de la molécula es dada dentro 
de la aproximación tight-binding con el parámetro de Hopping t,,, de modo que el 
hamiltoniano molécular es dado por H,, = — Y tm [(1) (4 +1| + |¿+ 1)(+]] siendo |¿) los 
estados cuánticos asociados a los orbitales tipo s de los sitios atómicos de la molécula. 
Del mismo modo la interacción entre el nanoporo y la molécula se asume tambien dentro 
del modelo del tight-binding, de modo que que cada sitio de la molécula interactúa con 
el primer sito mas próximo del nanoporo por medio del parámetro de hopping tip dado 


por 








RA ano + (5.12) 
55 amb. Dra) Tay | 
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siendo hb = 47 yd y la distancia de cada sitio de la molécula en la dirección paralela 
x y perpendicular y a la coriente eléctrica, hacia el sitio atómico mas próximo del borde 


que define el nanoporo. 


Los resultados obtenidos para el coeficiente de transmisión del nanoporo con la 
presencia de una molécula lineal de cuatro átomos ubicada en el centro del nanoporo 
como se muestra en la figura [53] es representada en la figura donde se considera 
que cada sitio atómico de la molécula interactúa con sus vecinos próximos de la propia 
molécula con el parámetro de Hopping t,, = 1.2t, las curvas de color negro representan 
el coeficiente de transmisión del sistema con nanoporo y la presencia de la molécula y 
la curva en color azul es el coeficiente de transmisión del sistema con nanoporo pero si 
la presencia de la molécula. La presencia de la molécula dentro del nanoporo produce 
resonancias y antiresonancias en el coeficiente de transmisión, que ocurren en cuatro 
valores de energía E = +2.05 eV ;+5.34 eV, estas resonancias y antiresonancias (ver 


inset de la Fig. corresponden al acoplamiento entre los sitios de la molécula. 





0 
Energia (eV) 


Figura 54: Coeficiente de transmisión de nanoporo rectangular de grafeno bajo la 


presencia de una molécula lineal de cuatro átomos de la figura 


Debido a que las resonancias y antiresonancias que se observan en los coeficientes 


de transmisión son pequeñas en las figuras siguientes representamos las variaciones del 
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coeficiente de transmisión AT' en lugar del coeficiente de transmisión 7 (E), AT' mide 
la diferencia del coeficiente de transmisión del sistema con nanoporo y el coeficiente de 
transmisión del sistema con nanoporo y la presencia de la molécula, de modo que las 
resonancias se representan por picos negativos en AT' y las antiresonancias por picos 
positivos en AT. Para las figuras siguientes consideramos un nanoporo de longitud 
D = 7y3a y ancho W = 7a, cuyo coeficiente de transmisión es representado en la 
figura |[52(e). Los resultados mostrados en la figura [55| representan AT' en función de 
energía del sistema con nanoporo y la presencia en el centro del nanoporo de moléculas 
lineales compuesta por distinto número de átomos. Los átomos que constituyen a la 
molécula se consideran idénticos e interactúan entre sí con un parámetro de Hopping 
t,, =t. 













































1 [ T T T T T ] T T T ] 1 
0.8- (a) 70.8 
06F 50.6 

a 0.4 7] 0.4 3) 
E 021 30.2 E 
0 
-0.2 
-0.4 1 1 | 1 | 1 | 1 | 1 -0,4 
Energia (eV) Energia (eV) 

ll y 7 1 
0.8- 70.8 
0.6 E 50.6 

a 0.47 30.4 ÉS) 
E 021 30.2 E 

0 0 
0.2 -0,2 
-0,4 -0.4 

-5 0 3 -5 0 5 
Energia (eV) Energia (eV) 


Figura 55: Variación del coeficiente de transmisión AT de nanoporos rectangulares de 
grafeno bajo la presencia de moléculas lineales formadas por diferente número de átomos 


identicos (a) 3 átomos, (b) 4 átomos, (c) 5 átomos y (d) 6 átomos. 


Los resultados obtenidos para AT' cuando la molécula en el centro del nanoporo es 
compuesta por tres átomos colineales se muestran en la figura B5la), donde se observan 
tres picos para AT en las energías E = 0; +3.9 eV; para una molécula lineal compuesta 
por cuatro átomos idénticos, la figura [55[b) muestra cuatro picos en AT' en las energías 
E = +1.8;+4.5 eV; si la molécula se compone de cinco átomos existen cinco picos en 
AT en las energías E = 0.0; +2.8; +4.8 eV como muestra la figura [55[c); y finalmente 
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si la molécula se compone de seis átomos AT' presenta seis picos, que ocurren en las 
energias E = +0.98; +2.78;+4.0 eV mostrados en la figura 5d). Estos resultados 
indican que la distribución espectral de la energía de acoplamiento de un número de 
átomos que componen la molécula lineal generan exactamente el mismo número de 
picos en el perfil de AT', asímismo la intensidad de los picos en AT' son mas intensos 


para determinados valores de energía. 
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Figura 56: Variación del coeficiente de transmisión AT' de nanoporo rectangular de grafeno 
bajo la presencia de una molécula lineal de cuatro átomos de la figura [53] con distintos 


parámetros de Hopping tm (a) tm = 0.6€t, (b) tm = 0.8€, (c) tm = 1.2t y (d) tm = 1.42. 


Con la finalidad de emular el efecto de moléculas formadas por cuatro átomos 
colineales de diferentes tipos, calculamos los perfiles de AT' considerando que cada sitio 
atómico de la molécula interactúa entre sí con parámetros de Hopping t,, variable. Los 
resultados se muestran en la figura [56] donde la figura |56[a) corresponde al parámetro 
de Hopping entre sitios de la molécula dados por t,, = 0.6t, resultando que los picos 
en AT ocurren en las energías E = +0.7; +1.8 eV; la figura [56[(b) corresponde a una 
molécula lineal con cuatro átomos que interactúan entre sí con el parámetro de Hopping 
tin = 0.8t encontrándose que los picos en AT' ocurren en las energías E = +1.0;,2.7 eV; 
si el parámetro de Hopping entre sitios de la molécula es t,, = 1.2t los picos en AT 


ocurren en las energías E = +2.0;+5.4 eV como muestra la figura [56[c), finalmente la 
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figura B6[d) corresponde al caso cuando los átomos de la molécula lineal interactúan 
entre sí con el parámetro de Hopping t,, = 1.4t apareciendo los picos en AT' en las 
energías E = +2.5;+6.4 eV. De acuerdo a estos resultados inferimos que el nanoporo 
permite leer moléculas con estructuras similares pero compuesta por átomos distintos, 


dado que los picos que aparecen en AT' ocurren para valores de energía totalmente 
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Figura 57: Variación del coeficiente de transmisión AT' de nanoporo rectangular de grafeno 


con presencia de moléculas de distintos tipos mostrados en la figura [49] respectivamente. 


Por otro lado con el fin de emular el efecto producido en el trasporte de carga 
del nanoporo debido a la presencia de moléculas diferentes a las moléculas lineales 
consideramos en el nanoporo de longitud D = 743a y ancho W = 7a moléculas 
de las configuraciones mostradas en la figura teniendo en cuenta que cada sitio 
atómico de la molécula interactúan entre sí con el parámetro de Hopping tm = t, 
los espectros correspondientes para AT' se muestran en la figura en los cuales 
se observa que los picos en AT' ocurren para diferentes valores de energía, los picos 
positivos representan antiresonancias, debido a que AT' es dado como la diferencia 
entre T (E) del sistema con nanoporo y el T (E) del sistema con nanoporo y molécula, 
las resonancias se reprsentan por picos negativos de AT', de acuerdo a los resultados 
mostrados las antiresonancias predominan sobre las resonancias en todos los casos, 


hecho que implica que la presencia de las moléculas dentro del nanoporo producifa 
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en general un incremento de la resistencia eléctrica y consecuentemente provoca la 
disminución de la corriente eléctrica del nanoporo, así mismo debido a las formas en 
las cuales fueron tomados los modelos de molécula el número de picos no coincide con el 
número de átomos de constituyen la molécula, pero los resultados son particularmente 


característicos del tipo de molécula considerada. 


Conclusiones 


En un nanoporo rectangular construido en una red finita de grafeno conectado a 
contacto ideales el transporte de carga depende del ancho del nanoporo, afirmación que 
se sustenta en las condiciones de transporte balístico considerado y que de acuerdo a los 
resultados el coeficiente de transmisión disminuye a medida que el ancho del nanoporo 


aumenta. 


Las moléculas presentes en el centro del nanoporo rectangular producen picos en 
el espectro de AT' que son característicos para cada tipo de molécula, dichos picos 
dependen del tipo de átomo, y el número de átomos que forman la molécula, asimismo 
depende de la forma estructural en el que cada uno de los átomos se distribuyen dentro 
de la molécula. En particular para moléculas lineales el número de picos de AT' coincide 
exactamente con el número de átomos constituyentes de la molécula, y varía con el 
tipo de átomos de forman la molécula, esta afirmación no se verifica en moléculas no 
lineales, pero se encuentra, que los espectros de AT' son muy particulares en cada 
caso. Las paticularidades que presentan los espectros de AT' abren la posibilidad de 
utilizar los nanoporos de grafeno como una técnica espectroscópica para el análisis de 
materiales, para tal fin se requiere construir una base de datos completa de cambios 
en la conductancia eléctrica (coeficiente de transmisión) producida por las diferentes 
moléculas que existen en la naturaleza y así mismo desarrollar tecnicas experimentales 


que permitan medir los pequeños cambios en el transporte de carga. 
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Capítulo 6 


Conclusiones 


En el punto cuántico con 4 contactos la fórmula de Landauer-Bitiker indica que 
la resistencia Hall Ry es asimétrico respecto a Ep = 0 a diferencia de la resistencia 
longitudinal R, que es simétrica respecto a Ep = 0. Para campos magnéticos débiles 
(a: pequeño) la formación de los niveles de Landau depende del tamaño del punto 
cuántico, observándose estos niveles con mayor definición cuando las dimensiones del 
punto cuántico son grandes. Asimismo el valor absoluto de la resistencia Hall |Ry| 
alcanza picos muy altos para determinados valores de energía en las cuales el campo 


magnético produce transisiones de fase cuánticas. 


El transporte electrónico de los nanoribbons depende del ancho de los nanoribbons 
y las condiciones de borde zigzag, armchair y/o bearded, ocurriendo que paracualquier 
ancho W los nanoribbons zigzag tienen comportamiento conductor, los nanoribbons 
armchair tienen comportamiento semi-metálico si el número de dímeros N de la cadena 
unidad verifica que M = (N + 1)/3 es un número entero positivo, caso contrario 
los nanoribbons armchair tienen comportamiento semiconductor. Los nanoribbons 
de borde bearded tienen comportamiento conductor, y presentan perfiles de 7 (E) 
idénticas a los nanoribbons zigzag. Los nanoribbons de borde mixto zigzag - bearded 
tienen comportamiento semiconductor. En los casos en que el comportamiento de los 
nanoribbons es semiconductor el gap de energía entre la banda de conducción y valencia 
tiene relación inversa con el ancho del nanoribbons, es decir, el gap de energía disminuye 


a medida que el ancho del nanoribbons incrementa. 


Un campo magnético uniforme que actúa perpendicularmente sobre un nanoribbons, 


introduce en los parámetros de Hopping un cambio de fase dado por la Sustitución de 
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Pierls, que produce la aparición de niveles de Landau que se observan cuando los 
campos magnéticos aplicados son débiles, estos niveles de Landau se transforman en 
subbandas de Landau cuando el campo magnético aplicado es extremadamente alto y 


el ancho del nanoribbons es mucho mayor al radio ciclotrónico. 


En un nanoporo rectangular construido en una red rectangular que emula al GaAs 
o en el grafeno el transporte de carga depende del ancho del nanoporo, afirmación 
que se sustenta en las condiciones de transporte balístico considerado y los resultados 
mostrados en las figuras y que indican que los coeficientes de transmisión 


disminuyen a medida que el ancho del nanoporo aumenta. 


La presencia de moléculas en el centro del nanoporo rectangular construida en na red 
rectangular o el grafeno producen picos en el espectro de AT' que son característicos 
para cada tipo de molécula, dichos picos dependen del tipo de átomo, y el número 
de átomos que forman la molécula. Asimismo los valores de energía en las cuales 
se presentan los picos de AT' dependen de la forma estructural en el que cada uno 
de los átomos se distribuyen dentro de la molécula. En particular para moléculas 
lineales el número de picos de AT' coincide exactamente con el número de átomos 
constituyentes de la molécula, y varía con el tipo de átomos que constituyen la molécula, 
esta afirmación no se verifica en moléculas no lineales, pero se encuentra, que los 
espectros de AT' son muy particulares en cada caso. Las paticularidades que presentan 
los espectros de AT' abren la posibilidad de utilizar los nanoporos de grafeno como una 
técnica espectroscópica para el análisis de materiales, para tal fin se requiere construir 
una base de datos completa de cambios en la conductancia eléctrica (coeficiente de 
transmisión) producida por las diferentes moléculas que existen en la naturaleza y así 
mismo desarrollar tecnicas experimentales que permitan medir los pequeños cambios 


en el transporte de carga. 


De la observación de las resonancias y antiresonancias producidas por la presencia 
de moléculas en el centro de nanoporos rectangulares en red cuadrada y en el grafeno, 
concluimos que en nanoporos construidos en grafeno con la presencia de moléculas las 
antiresonancias son predominantes respecto a las resonancias, en cambio en nanoporos 
construidos en redes cuadradas la predominancia de las resonancias o antiresonancias 
depende de la orientación de la molécula respecto al sentido de la corriente del 


nanoporo. 


Los métodos numéricos de decimación y recursivo de convergencia rápida son 


eficientes para evaluar el transporte electrónico utilizando la formula de Landauer, 


102 


siendo el método recursivo de convergencia rápida mas eficiente que el método de 
decimación para determinar las funciones de Green de los contactos, pero para el 
análisis del nanodispositivo en sí es más recomendable el método de decimación por 
ser aplicable incluso cuando las cadenas verticales consideradas tienen estructuras 
diferentes, casos en las cuales el método recursivo de convergencia rápida no es 
aplicable, es decir una condición necesaria para la utilización del método recursivo 
de convergencia rápida es que las cadenas verticales que constituyen el nanodispositivo 


sean exactamente iguales. 
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Apéndice A 


Relaciones importantes 


En la sección 2.2 se utilizan matrices hamiltonianas H que permiten obterner la 
estructura de bandas, aqui derivamos explícitamente las expresiones de los elementos 


de matriz. 


A.l. Elementos de matriz para el hamiltoniano del 


grafeno 


En esta sección calculamos los elementos de la matriz hamiltoniana H dada en 


la ecuación (2.11) en la base (|), |1)), dadas en (2.7), (2.8) y de acuerdo al 
hamiltoniano tight-binding (2.3), es decir el elemento de matriz (4*|H |) se obtiene 


como: 


(Hat) = Ne REA | PI + 107) 61) | Y Rm óA) 
l (1,5) m 
= Verlo + (tr o | AA oa) 
l (4,5) m 


y is 3 ME y 95 ó?| y er 16%,) 
l (1,5) m 
= AYMARA PO) =0 
l (a) m 
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(01H IE) =0 (A.1) 


de manera similar se encuentra 


(15 1 |") = 0. (A.2) 


Calculando ahora el elemento de matriz (| H|4P) tenemos 


(oe EE) 


SAR All ES (16681 +107)(68 1) | ER O) 
l 


(1,5) me 


Ne ER? [Y (01070621 + (01 167) DS ) 
l (6,5) 


A 20 ón] A 
=p a era, 
AAN Y et 
—t y ¿ER y y AR 
= —t y ener y E 
sq y y ¿RPRE) 

la) 


donde RP — Rf representa la posición relativa de los átomos de la red B respecto a 


las posiciones de los átomos A. 


Para un solo átomo A tenemos que las posiciones relativas de los primeros vecinos 
son RP — Rf = 61,02, 03, de modo que 


WwéHIP) = EI == [est y ¿by ab) 


(eE [wE) =—tf(K) (A.3) 
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con F(k) = cikó + cikd + cikó3 
Por otro lado 
(pH) 
A eras 
(a 
+ yo y ¿k(RÍ-RP) 
(a 


(0 Hao) 


para un solo átomo B las posiciones relativas de los primeros vecinos son Rf — R* = 


—01, —02, —03, entonces 


(eE (wí) = as ik (RÍZRL) — _ = + (0 ME) y ¿ED 4 gen) 


(e/ Hi) =—t de só e = —tf*(K) (4.4) 


Sustituyendo los resultados A.1), (4.2), (A.3 y en la ecuación (2.11), la 


expresión matricial de H es 
k 
H=-—t 2 FU) z (A.5) 
Fx) 0 


A.2. Relación de dispersión |f(k)| 


Evaluamos la relación de dispersión |f(k)| teniendo en cuenta que f(k) es dado 
como f(k) = ele% 4 gikóz 4 gikós es decir 


109P = PU) o 
IF()]2 = (e7ieó q grikeós 4 grileds (gileór y gileds 1 gileds) 








3+e ik-(Í1 te ik-(81 Be ik-(d2 De ik-(d2 te ik-(S3 De ik-(S3 82) 
= 34% (8,82) y ek 82) y pik (8183) y e k(81 33) y ei (82—33) y e 7ik:(82—33) 
3 + 2cos(k - (8, — 92)) + 2cos(k - (8, — 93)) + 2cos(k - (92 — d3)) 
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do 3. 1 3. 1 
a Pas Ñ 205) A (Pos E 205) = —V3a4 


d — Ó 


1! 
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| 
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Sa 
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| 
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>> 
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Sa 
3 
| 
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! 
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[5 
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>> 
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como cos(—0) = cos 9 paridad de la función cos 9 





f(x) 1? = 3 + 2cos (V3k.a) + 2c0s (Pra ha) + 2c0s (Pa 5 ha) 


usando cos (p + 0) = cos Hcos Ó + sin psin O) 


COS (Pa + 3 ya = COS 7 Ka cos (Gta) — sin 7 Ka sin (Gta) 
Cos (hs + 3 ya = Cos 7 Ka COS (Gta) - sin 5 k, a | sin (Sta) 


entonces 








|F(k)1? =3 +2 cos (V3k,a) + 4 cos (Pra) COS (Gta) (A.6) 


si cos? ($) = 3 (1 +c0s 20) entonces cos 29 = 2cos? ($) — 1 


cos (Vta) =D608* (Pra) =-1 (A.7) 
2 
ó a v3 3 
|F(k)// =1 +4c0s (Pra) + 4cos (Pra) COS (Eta) (A.8) 


de donde 





|F(k)| = 1 +4cos (Pra) COS (Gt) + 4 cos? (Pra) (A.9) 
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